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ABSTRACT 

We present the first systematic investigation of spectral properties of 17 Type Ic Supernovae (SNe Ic), 
10 broad-lined SNe Ic (SNe Ic-bl) without observed Gamma-Ray Bursts (GRBs) and 11 SNe Ic-bl with 
GRBs (SN-GRBs) as a function of time in order to probe their explosion conditions and progenitors. 
Using a number of novel methods, we analyze a total of 407 spectra, which were drawn from publishe d 
spectra of individual SNe as well as from the densely time-sampled spectra of Modjaz et al. (2014). 
In order to quantify the diversity of the SN spectra as a function of SN subtype, we construct average 
spectra of SNe Ic, SNe Ic-bl without GRBs and SNe Ic-bl with GRBs. We find that SN 19941 is 
not a typical SN Ic, in contrast to common belief, while the spectra of SN 1998bw/GRB 980425 are 
representative of mean spectra of SNe Ic-bl with GRBs. We measure the ejecta absorption and width 
velocities using a new method described here and find that SNe Ic-bl with GRBs, on average, have 
quantifiably higher absorption velocities, as well as broader line widths than SNe without observed 
GRBs. In addition, we search for correlations between SN-GRB spectral properties and the energies 
of their accompanying GRBs. Finally, we show that the absence of clear He lines in optical spectra 
of SNe Ic-bl, and in particular of SN-GRBs, is not due to them being too smeared out due to the 
high velocities present in the ejecta. This implies that the progenitor stars of SN-GRBs are probably 
He-free, in addition to being H-free, which puts strong constraints on the stellar evolutionary paths 
needed to produce such SN-GRB progenitors at the observed low metallicities. 

Subject headings: supernovae: general - Gamma-ray bursts: general - supernovae: individual: SN 
1992ar, SN 19941, SN 2004aw, SN 2007gr, 2011bm, PTE 12gzk, SN 1998bw, SN 
2003dh, SN 20031w, SN 2006aj, SN 2009bb, SN2009nz, SN 2010bh, SN 2010ma, 
SN2012ap, SN 2012bz, SN 2013cq, SN 2013dx, SN 2013ez — gamma-ray burst: 
individual: GRB 980425, GRB 030329, GRB 031203, GRB 060218, GRB 100316D, 
101219B, GRB 120422A, GRB 130427A, GRB 130702A, GRB 130215A 


1. INTRODUCTION 

Stripped supernovae (SNe) and long-duration Gamma- 
Ray Bursts (long GRBs) are amongst nature’s most pow¬ 
erful explosions from massive stars. They energize and 
enrich their host galaxies, and, like beacons, they are 
visible over large cosmological distances. However, their 
progenitor systems remain to be identified and the details 
of the explosion mechanis m ar e unknown (se e reviews by 
Woosley & Bloom||2006| and |Smartt||2009|). Stripped- 
enveiope SNe (i . e, oJNe of Types iTR lb, and Ic , e.g., 
Filippen]^|1997 Modjaz et al.||2014 |, in short “stripped 
bNe," are core-collapse events whose massive progeni¬ 
tors have been removed of progressively larger amounts 
of their outermost H and He envelopes. 

Specifically, broad-lined SNe Ic (SNe Ic-bl) are SNe 
Ic whose line widths approach 15,000—30,000 km s~^ 
around maximum light and whose optical spectra show 
no obvious presence of H and He. Nevertheless, the ex¬ 
act distinction between a normal SN Ic and a SN Ic-bl 
is still debated, as well as whether there is hidden He in 
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SN Ic-bl spectra. 

For the last fifteen years, the exciting connection be¬ 
tween long GRBs and stripped SNe has been observation- 
ally well established, in particular with SNe Ic-bl, as the 
only type of SNe observed accompanying long GRBs for 
both low-luminosity, nearby, an d high-luminosity, cosmo¬ 
logical GRBs (for reviews, see Woosley & Bloom 2006 


[Hjorth &: Bloom|2011[ |Modjaz|2011| and tor the most re^ 

cent SN accompanying a high-luminosity GRB,prn et al.| 


2013). However, the existence of many more SNe Ic- 
bl without GRBs (most convincin gly SN Ic-bl 2002ap 
without a relativistic engine, e.g., Berger et al.||2002| 
raises the question of what distinguishes SJN-GRB pro- 
genitors from those of ordinary SN Ic-bl without GRBs. 
Viewing angle effects are probably not the only reason 
why those SNe Ic-bl did not s how an accompanying GRB 
(Soderberg et al.||2006 2010). One promising line of re- 
search is to investigate what sets apart SNe with GRBs 
from those without GRBs in order to elucidate the con¬ 
ditions and progenitors of these two types of explosions. 
While there are numerous possible avenues (for recen t 
reviews, see e.g., Hjorth & Bloom 2011 Modjaz 2011), 
here we concentrate on comparing the optical spectra of 
SNe Ic-bl with GRBs to those without GRBs, as well as 
to those of normal SNe Ic. Optical, early-time, i.e., pho- 
tospheric phase, spectra are used for the classification 
of different explosions and probe the bulk of the ejected 










































2 


stellar material, in particular the outermost layers. Prior 
work involving syntheti c models based on Monte Carlo 
radiative transfer codes (Mazzali et al. ]M^ [Walker et al. 
'2014 and references therein], while important, has in¬ 
cluded only two normal SNe Ic (SN 19941, SN 2004aw) 
and only four SN Ic-bl without GRBs, thus not yet pro¬ 
viding a statistically significant sample. In particular, in 


other comparisons i n the literature (e.g.,[Pian et al. 2006 
Valenti et al.||2008 Taubenberger et al. |2006 ), SIN 1994. 


is assumed to represent the class of SNe Ic, while there is 
strong evidence that SN 19941 is not a “typical” SN Ic, 
neither pho tometrically (Richardson et al. 2006 |Drout 
et al. 20111, nor spectroscopically, as we will show be¬ 
low. 

The time is ripe to undertake this detailed compari¬ 
son, given that from the liter ature and as part o f the 
CfA Stripp ed SN data release (|Modjaz et al.|2014| here¬ 
after M14), we can now include a statistically significant 
set of 23 Si4e Ic-bl (with and without GRBs) and 17 SNe 
Ic that have accompanying photometry, so that we can 
assign phases to the spectra. Now, for the first time, 
we have a large enough sample to conduct a systematic, 
data-driven study on quantifying the similarity and be¬ 
havior of SNe Ic-bl (especially those accompanied with 
GRBs) compared to those with SNe Ic. Here, we concen¬ 
trate on the SN-GRB connection and use the same novel 


methodology we develop in Liu et al. (20151, in which we 
concentrate on SNe lib, Ib and Ic. 

The format is as follows: in Section [2 we summarize 
the sample of SNe used in this paper, "raich we analyze 
in the remainder of the paper using two complimentary 
methods. In Sectionwe quantify the spectral diversity 
of stripped-envelope sNe by analyzing and comparing 
mean spectra, as well as corresponding standard devi¬ 
ation and maximum deviation spectra, which we con¬ 
struct for SNe Ic, SNe Ic-bl, and SNe Ic-bl connected 
with GRBs. This method takes into account all features 
that compose the spectrum. We also use mean spectra 
to explore the question of whether SNe Ic-bl may have 
smeared out He lines hidden in their spectra - one of the 
outstanding questions in understanding the progenitors 
of SN-GRBs as it is hard to produce a star that is both H- 
and He-free at the low metallicities at which SN-GRBs 
are observed. In Section]^ we concentrate on the velocity 
evolution of the Fe II A5169 line and conduct a compari¬ 
son of its behavior (both in terms of its Doppler velocity 
and its width) amongst the SN types, for which we devel¬ 
oped a new method for measuring blended lines in spec¬ 
tra of SNe Ic-bl, described in the Appendix. In Sectionj^ 
we discuss our results and their implications for progeni¬ 
tor models of SNe Ic-bl, specifically those with and with¬ 
out observed GRBs, and conclude with Section]^ 

2. SN SPECTRAL SAMPLE AND OTHER RELEVANT 
PROPERTIES 

We list our SN spectral sample and their references in 
Table for normal SNe Ic, and in Table for SNe Ic- 
bl witn and without observed G RBs. We only incl ude 
SNe with secure spectral IDs (see Modjaz et al.|2014 and 
below for more information), whose spectra have been 
published before February 2015 and were made available 
to us by the authors, and whose date of maximum light 
has been measured. With those criteria, we have in¬ 
cluded 17 SNe Ic and a total of 23 SNe Ic-bl, both with 


and without observed GRBs. 

Amongst the latter group, we have included 11 SNe Ic- 
bl with GRBs, sometimes referred to as SN-GRBtH with 
published spectra, and denote them as “SNe Ic-bl with 
GRBs” in plots. For the SN-GRBs, we further distin¬ 
guish SNe Ic-bl accompanied by low-luminosity GRBs 
(LLGRBs) as opposed to those with high-luminosity 
GRBs (HLGRBs), where the luminosity refers to the 
isotropic gamma-ray luminosity. LLGRBs, in addition 
to their low luminosities (lO'^® — 10^® erg s“^), ap¬ 
pear to be more isotropic and have relatively soft high- 
energy spectra compared to HLGRBs (e.g., Gampana 
|et al.| |2006[ [Hjorth |2013[ |Nakar| |2015| and references 
thereinb There have been suggestions that the high- 
energy emission of the GRBs in these two groups may 


come from different mechanisms (e.g., Waxman et al. 
2007 1 [Bromberg et al.|[2(311[ [Nakar fc bari| ^12[ |1N akar 
20151). We note that when corrected for beaming an- 


gle, the intri nsic energies of LLGR Bs and HLGRBs are 
similar (e.g., Mazzali et al. 2014), and that there are 
also some GRBs that appear to be transition objects 
(e.g., SN 2012bz/GRB 120422A, iSchulze et al.l[20T4) 
with intermediate luminosities, which leads us to col- 
lect the gamm a-ray luminosities of all the SN-GRBs 
in this sample (Mazzali et al. 2014). In Table [2 we 
mark the following LLGRBs: SN 1998bw/GRB 9^425, 
SN 20031w/GRB 031203, SN 2006aj/GRB 060218, SN 
2010bh/GRB 100316D. We furthermore analyze their SN 
spectra to see if there are differences between the spec¬ 
tra— features of the SNe connected with these two kinds 
of GRBs, and whether there are any correlations between 
SN spectral properties and GRB en ergies (all GRB en- 
ergies taken from the compilation in Mazzali et al.[[20T4 
see references therein). 


et al. 


In addition, we included SNe I c-bl 2009bb (Soderberg 


et al. 


Mlllsavlj evic et al. 


20101 iPignata et al.[[2011| and 2012ap (jiVIarKutti 

mf[Mi" .. . " • ^ 


015[) as “engine-driven” 
SNe since they were suggested to be relativistic SNe 
based on their large and early radio emission, similar to 
nearby SN-GRBs. even if no gamma rays were d etected 
([Soderberg et al.[|2010| [Ghakraborti et al.[[2015[). Thus, 
in our analysis below, we include SNe 2U09bb and 2012ap 
in our sample of SN-GRBs. 

We also included SNe Ic-bl without detected gamma- 
ray emission and designate them as “SNe Ic-bl with¬ 
out observed GRB” in the subsequent analysis, while 
referring to them simply as “SNe Ic-bl” in the plots 
(as opposed to “SNe Ic-bl -I- GRBs”). The presence 
of an off-axis or weak jet can be excluded for some, 
but not all, SNe Ic-bl without observed GRBs in our 
sample. For example, while for SN Ic-bl 2002ap, the 
radio observations were deep enough t o exclude the 
presence of a jet, be it on- or off-axis (Berger et al. 


We were not able to include one SN-GRB, since the 
authors did not provide u s with the i r pub lished spectrum: 
SN 2008hw/GRB 081007 (|Jin et ah] |2013|l. In addition, 
there are other SN-GRBs with spectra that have been re¬ 
ported in circulars, but are not included here as their 
spectra have not been pu blished yet: GRB 111211 A ij de 
Ugarte Post igo et al.[[2012[l. SN 2012eb/ GRB 120714B l[Klo^ 
let al.| [20121, SN 2Ul3iu/G ' ' 

■PTFR 

last stages of prep aring this manuscript, "the submitted paper by 
|Gano et al.| (|2015| with spectra appeared on arXiv after February 
2015, and thus, is not included here. 


.II 2 OI 2 L SN 2uV31u/ GRB 130831A (iK lose et all 
rabtu/GRB 140606B | |Perley et al.||20l4[ l, tor which in tne 
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TABLE 1 
Sample of SNe Ic 


SN name 

redshift 2 

Phases of spectra with respect to maximum light^ 

References'^ 

SN 1983Vt 

0.0054 

-15, -13, -10, -2, 6, 7, 8, 9, 10, 11, 12, 32 

SNID(C97) 

SN 1990B‘' 

0.0075 

4, 5, 6, 6, 7, 9, 10, 15, 22, 28, 28, 31, 39, 47, 58, 65, 72, 88, 90-l-(2) 

SNID(B90, MOl, CfA) 

SN 1992ar'' 

0.1450 

3 

COO 

SN 19941 

0.0015 

-6, -6, -5, -4, -3, -3, -2, -1, 1, 2, 2, 3, 3, 6, 7, 11, 21, 22, 23, 24, 

25, 26, 28, 30, 32, 33, 34, 34, 36, 37, 38, 40, 55, 56, 68,90+(3) 

SNID(F95, C96, M99), M14 

SN 2004aw 

0.0158 

-5, -3, -2, -2, -1, 0, 1, 2, 2, 3, 3, 4, 5, 12, 18, 18, 19, 21, 23, 

25, 26, 30, 32, 46, 48, 90-|-(2) 

SNID(T06), M14 

SN 2004dn 

0.0126 

-8, 31, 34, 52, 54 

H08, M14 

SN 2004feT 

0.0179 

-7, -7, -6, -5, -1, 0, 1, 2, 3, 9, 22, 90-l-(l) 

M14 

SN 2004ge 

0.0155 

12 

M14 

SN 2004gt 

0.0055 

16, 18, 19, 22, 43, 48, 70, 81, 90+(5) 

M14 

SN 2005az 

0.0088 

-8, -7, -5, -3, -1, 1, 4, 17, 18, 21, 23, 25, 27, 29, 30, 31, 32, 

48, 56, 61, 64, 86 

M14 

SN 2005kl 

0.0035 

-4, 39, 70, 90+(l) 

M14 

SN 2005mf 

0.0189 

-2, -1, 5, 8 

M14 

SN 2007cl 

0.0222 

-6, -6 

S12, M14 

SN 2007gr 

0.0017 

-9, -8, -6, -5, -4, -3, 0, 5, 7, 8, 9, 10, 13, 14, 15, 16, 17, 18, 

20, 25, 42, 45, 47, 49, 53, 53 

V08, M14 

SN 2011bmt 11 

0.0221 

-11, -3, -1, 16, 24, 30, 50, 61, 90+(5) 

V12 

PTF12gzkt t ttd 

0.0138 

-10, -5, -4, 0, 2, 3, 4, 6 

B12 

SN 2013dk 

0.0055 

-1 

E13 


® Phases are in the rest-frame with respect to P-band maximum light, either directly measured (as listed in |Modjaz et al. ]l(20^ or 
in the original papers or using the polynomial fit as described in|Bianco et al.|(|2014||), or transformed (see text tor details). l-'Uases 
are rounded to the nearest whole day. The number in a bracket is the number ot spectra with phases larger than 90 days after the 
date of maximum light - spectra which we include for completeness, but do not analyze here . 

^ References: SNI P = in SNIP release ver sion 5.0 via templates-2.0 by Blondin Sz Tonry ( 20071 . with the original references in 
parentheses; C97 = |Clocchiatti et al,' 

CfA S upernova Program 
C96 "" 

V12 = 

^ Note 


Clocchiatti et al, 

Valenti et ai. TO 


(|1997||; B9 0 =[Barbon et al.| ||1 990ll; MU l =[Matlieson et al 

jvi 99 — 


c publ ished he re; Ml4 =|Modlaz et al. jl2Uf4|l ; M9 
111996 I; H0 8 =|HarutvunWn et al .lJ^UU^; '1U6 =r" 
1^; l3l2 =|Beri-Ami et al.| l |2Ul2 i; liild ias-Iu 


= riaubenperger et a 


er ger et ^ 

.|||2U13l. 

their OPS 


h 


|2001| l ; CfA = spectra from th e 


iFilioDenko et al. 

1995 

= |ailverman et al. 

2TTr2 


.osa et al 

ttiat tlie dates ot maximum tor blNe iyyun and 11^92ar are uncertain since their observed light curves start after maximum, 
and their dates were estimated via fitting with the light curve of other well-observed Stripped SNe. 

^ While SNIP finds good matches of PTF 12gzk with SN Ic 2004aw, a SN Ic, w e note that it displays v ery high absorption velocities 
, similar to th ose of SN Ic-bl, as inferred from the optical spectra (e.g., |Ben-Ami et al.||2012|| and from radio observations 
3sh et al.|2013||. ^hus, while we call it a ”SN Ic” in the table, we are not including t' L t L'ZgzK wnen constructing mean spectra 
ectionj^for details). 


Horesh 
ot SJN Ic (see 

^ For SN 1983V, we determine the date of V-band maximum by applying the polynomial fit as described in 
the original data of C97, to be consistent with the date of maximum determination for the rest of the SNe iii our sample, 
of maximum is JDy = 2,445,681.4 ± 1.0, which is 6 days later than th e date obtained in C 97 with light curve template fitt ing. 


Bianco et al. 


(2014| 

Jut 


to 

ate 


Bianco et al. 


(20141, which is different from that in 


For SN 2 004fe. we use the date of V-band maximum as measured by 

ITT por 2011bm, we determine the date of V-band maximum by applying the polynomial fit as described in 
the original data of V12, and obtain JDy = 2,455,673.7di 2.1. 

11 tt Por PTF12gzk, we determine the date of V-band maximum by applying the polynomial fit as described in 
to the original data of B12, and obtain JDy = 2,456,150.1 di 0.7. 


Prout 


Bianco et al. 


120141 to 


Bianco et al. 


(20141 


20021, for SN 2010ay, Sanders et al. (2012b) cannot 
rule out a low-luminosity GRL5 like GRSOOohs from 
their radio and X-ray upper limits. In general in all 
SNe Ic-bl without observed GRBs, except in the most 
nearby SN 1997ef, 2002ap and SN 2003jd, the various 
X-ray and radio l imits, when t hey exist (for PTF lObzf: 
Corsi et al. 2011[ PTFlOvgv: Corsi et al. 2012), only 


rule out very energetic relativistic UKB explosions, such 
as SN 2003dh/GRB 030329, but not necessarily low- 
luminosity GRBs such as SN 1998bw/GRB 980425 or 
SN 2006aj/GRB 060218. With the caveat that the pres¬ 
ence of off-axis and low-luminosity GRBs cannot be ruled 
out for all SNe Ic-bl in our sample, we proceed in the 
rest of the paper to designate these as SNe Ic-bl without 


GRBs, but discuss this caveat in more detail later on. 

Unless found by the accompanying GRB trigger, the 
SNe in our sample were discovered by both galaxy- 


targe ted SN surveys (e.g., LOSS, Filippenko et al. 


2001) a nd galaxy-unta rgetted ones (e.g., 'T exas SN 


Seardl^ 

Quimby 

20061 

et al. 

|2008p. Here, we 


,2ic-bi=0.0232 for the SN Ic-bl without observed GRBs, 
and zsNGRB =0-1450 for the SN-GRBs. The fact that the 
different SN subtypes span somewhat different redshift 


http: //grad40. as .utexas . edu/~quimby/tss/index.html 
® http://WWW.sdss.org/supernova/aboutsprnova.html 
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TABLE 2 

Sample of SNe Ic-bl 


SN name 

Redshift 2 

Phases of spectra with respect to maximum light^ 

References'^ 

SN 1997ef 

0.0113 

-14, -12, -11, -10, -9, -6, -5, -5, -4, 7, 13, 14, 16, 17, 19, 
20, 22, 27, 38, 40, 44, 46, 48, 48, 75, 80, 86, 90-l-(l) 

100, M14 

SN 2002ap 

0.0022 

6, -5, -2, -1, 0, 1, 2, 3, 4, 5, 5, 6, 7, 7, 10, 12, 13, 

26, 27, 30, 31, 31, 90-l-(12) 

SNID(G02, F03), M14 

SN 2003jd 

0.0188 

-3, -2, -1, 0, 1, 19, 20, 21, 22, 23, 24, 27, 28, 34, 

46, 47, 48, 49, 49, 51, 54, 55, 60, 73, 90-l-(l) 

VOS, M14 

SN 2007D 

0.0232 

-8 

M14 

SN 2007bgt 

0.0346 

6, 7, 19, 27, 59 

YIO, M14 

SN 2007m 

0.0155 

-3, 0, 4, 5, 9, 13, 16, 29, 36, 59, 90-l-(l) 

S09, M14 

SN 2010ayt 

0.0671 

16, 26 

Sal2 

PTFlObzf (SN 2010ah)t 

0.0498 

2, 7 

Cll 

PTFlOqtst 

0.0907 

-1, 1, 10, 17, 20, 24, 90-1(1) 

W14 

PTFlOvgvt'^ 

0.0142 

-7, 4, 8, 37, 73 

C12 

SN 1998bw/GRB 980425* 

0.0085 

-11, -9, -8, -5, -4, -3, -1, 1, 2, 4, 7, 9, 10, 11, 17, 

20, 27, 43, 50, 62, 71, 90+(5) 

G98, POl 

SN 2003dh/GRB 030329 

0.1685 

-10, -9, -8, -7, -6, 6, 12, 14, 18, 30, 56 

S03, M03, K03, K04, D05 

SN 20031w/GRB 031203t* 

0.1006 

-2, 7 

M04 

SN 2006aj/GRB 060218* 

0.0335 

-7, -6, -6, -5, -5, -4, -3, -3, -2, -2, -1, -1, 0, 0, 

1, 2, 2, 3, 3, 4, 5, 6, 8, 9, 10, 90-1(1) 

M06, P06, M14 

SN 2009bb'i 

0.0099 

-5, -4, 1, 5, 15, 16, 20, 21, 24, 28, 42 

Pll 

SN 2009nz/GRB 091127® 

0.490 

~ 1 

Bll 

SN 2010bh/GRB 100316Dt* 

0.0593 

-4, -4, -4, -4, -3, -2, -1, 0, 0, 1, 1, 2, 4, 5, 6, 

8, 10, 12, 14, 15, 17, 24, 30, 33 

CIO, B12 

SN 2010ma/GRB 101219B 

0.55 

-13, -3, 10 

Sll 

SN 2012ap‘i 

0.0121 

-5, -3, -1, 1, 2, 4, 7, 8, 14, 25 

M15 

SN 2012bz/GRB 120422A 

0.2825 

1, 4, 9 

S14 

SN 2013cq/GRB 130427At 

0.3399 

-1 

X13 

SN 2013dx/GRB 130702A 

0.1450 

-5, -5, -3, -2, -1, 1, 1, 2, 5, 5, 6, 11, 13, 17, 20, 23 

D15 

SN 2013ez/GRB 130215Af 

> 0.597 

~0 

C14 


Phases ; 

after the date of ma ximum light • 
^ References: 100 = 


I day. 


days 


bvlBlo ndin &: Tonry 

"OF — 


Iwamoto et al. 
|2UU7||. with ttie' 


spectra with phases larger i 
sp ectra which we include for complete ness, but do not analyze here. 

(2000i; M14 = Modjaz et al. (2014l; SNI P = in SNIP release v ersion 5. 0 via templates-2. 0 
original references in pare ntneses;; G02 = Gal-Yam et al.| l|2002||: F03 = [Foley et al.| (|2003|[; 


,, _, . MU3 =1 

let al.lJ ^UUbll; iVlC 4 = IM alesani e 
(jiU == | u]iornock et al.N! 
[Schulze et al.| (|2U14||; Al 


W14 = 


VOa = I Vale nti et al .lil^Un^; YIO =|Y oung et al.l ||2Q10||; S09 = ISa hu et al. 
I^ T^tefll: (Ji2 ^ l U orsi It al.'^ . - ' ■ ■ ' . 

= [atanek et ah] ( 200^ 


]S2om 

[Mat ties 


latti eson e t al. 
et al.|l|2UU4||; M(J^ 
I2010fc: B12 =lBulano et ~ 


Walker et al.| ( |20l4|l. = 


_ .. _ = ^u et al.f(|2ULi|): Ui4= 

While [Uorsiet ak] l|2012[ | classity Fi'hlUvgv a SN 


(|2U U15ll; KU;j ^^Kawabat a et 
1 =|Modiaz et akl(|:^uuh|l; J^L 


4odia z et al.] ( 2UUb||; J-'Ub = 
2UI2I1: ail == bparre et al.' 
ano et ahl ||2Ui4|l. iJ15 =f 
ic based on its measureci 


Ulj|); Sai:^ = Sail ers et a i.j (|2ui:^at: uii = Gofsi] 
Galarna e t al.|l|ilj^ a|): FUl ^iPatat et a l.| (2001 : bUd 
in l IMi^ ; KDT^TKo sugi et al.[ l[M)4|rDbF^ P^ 
b == Jr'ian et al.l (|2 UUbl): J-"ii =[Fignata et al. l (|2uil|l; 
al. J2Ull|): Mio =|Milisavljevic et al.| l|2Ui^; ai4 = 


hllia et al.| 

low bi 11 Ab^bb absorption velocities, its optical 
with SNIP finding good matches with only SNe Ic-bl. Thus we 


spectra show broad lines amd are clearly those of a SN Ic-bl, 

re-classify this object as a SN Ic-bl, following our classification approach outlined in M14. 

^ We include SNe 2009bb and 2012ap as ’’engine-driven” SNe Ic-bl, since they have been suggested to be relativistic SNe based 
on their large radio-emission, sim ilar to, albeit weaker than, nearby SN-GRBs, even if no gamma-rays were detected ( |Soderberg| 


|et al.|2010| iMargutti et al.|2014||. _ 

Lor blN 2UUynz/GH15 Uyili7,|llerger et al.|||2011|) did not provide us their published spectrum, so we had to digitize their Figure 
1, whic h shows their smoot hed spectrum. Thus, some of our analysis which requires the raw spectrum, could not be performed. 

^ While Cano et al.[ (20141 suggest that SN 2013ez/GRB 130215A is a SN Ic, SNIP finds good matches with both SNe Ic-bl 


Cano et al.[(2014i suggest that SN 2013ez/GRB 130215A is a SN Ic, SNIP finds good matches with both SNe Ic-bl 
(including SiNe Ic-bl 2UU7ce and 2005nb ) and with SNe Ic. Thus, we call it a ”SN Ic/Ic-bl”, since SNIP does not find an 
unambiguous type, and we do not include it when constructing mean spectra of SNe Ic-bl. Furthermore, the redshift of this 
SN/GRB is uncertain, as th ere is only a lower limit f rom m etal absorption lines superimposed on the GRB afterglow spectrum 
as mentioned in the GCN by |Cucchiara &; Fumag alli| (|20 13||. 

^ For consistency, we used the transformations of Bianco et al.|||2014|l to convert the dates of maximum in the B- or R-band into 
that in the Y-band, as no Y-band data have been provided. 

These SN-GRB harbored low-luminosity GRBs. 
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ranges (with SN-GRBs at the highest redshifts, given 
their rarity) should not influence our results, as all SN 
spectra are de-redshifted into the rest frame, and all spec¬ 
tral comparisons are made for the spectral features that 
were observed in all SN spectra at all redshifts. In ad¬ 
dition, any evolutionary effect is expected to be insignif¬ 
icant within such a small redshift range, for which e.g. 
the metallici ty content of the univ erse does not change 
significantly ( Tremonti et al.|[2004 |. 

We note that we did not include a new and very rare 
kind of SNe Ic - the Superluminous SNe Ic, which are up 
to 100 times more luminous than normal core-collapse 
SNe and for which the powering source is high ly debated 
(e. g.,|Quimby et al .|2011||Chomiuk et al.|20111 see review 
by Gal- Yam| 2012). 'i'hey are not included here, since it 


is not yet clear how they are physically related to the 
“normal” SNe Ic and SN e Ic-bl in our sample (see dis¬ 


cussion m 


Gal-Yam 2012). There is a recent claim of a 


spectrum that looks like that of a super luminous SN Ic 
has been observed in connection with an ultra-long GRB 
(Greiner et al.|2015), though the spectrum is very noisy. 


2.1. SN spectra and their spectral types 

The spectra were drawn from p ublic archives (SUS- 
PECT, now subsumed by WISeREh ^[Yaron fc Gal-Yam 
2012), as well as the CfA Supernova Archiv^ with a 
large fraction of SNe Ic published in M14 or we directly 
requested them from the authors in the listed references. 
In order to concentrate on the SN features alone, we 
removed superimposed emission lines due to host HII re¬ 
gions. We also removed any telluric absorption in the 
non-Berkeley and non-CfA spectra if it was relatively 
narrow and a simple interpolation could reasonably re¬ 
move it (for more detai ls on pre-processing the SN spec¬ 
tra, see Liu et al.]|2015[ )p] 

We also include all relevant SNe that had b een initially 
included in t he SNID release version 5.C[^by Blondin fc 
Tonryl (20071, though in some cases with a nio difled SN 
identihcation (see extensive discussion in M14|. We also 
performed our own spectral classification on a fe w of the 
SNe using our expanded set of S NID templates (Modjaz, 
et al.|[2M4 |Liu fc Modjaz| 20I4) and found two cases m 
which the iL) was ditterent trorn that announced by the 
authors. 


Eor PTE lOvgv (Corsi et al. 2012), we find a better 
match with SNe Ic-bl than with SNe Ic, with SNID 
finding very good matches with a number of S Ne Ic- 
bl (SNe 199 7ef, 2006aj, 20031w, 2007ru). While [C^ 
et al. (|2012 ) classify it a SN Ic based on its measured 
Si II Ab355 velocities (which we also measure to be low, 
see below), its optical spectra are clearly those of a SNe 
Ic-bl, as also mentioned in the initial discovery ATel by 
PTE (Gorsi et al.||2010). Thus, following our SN classifi¬ 
cation approacFoutlined in M14, we reclassif y PTFlOvgv 
as a S N Ic-bl. For SN 2013ez/GRB 130215A,|Gano et ah' 
(2014) suggest that is a SN Ic based on their spectrum 
taken 16.1 days (rest-frame) after explosion, which cor- 


® http://www.weizmann.ac.il/astrophysics/wiserep/ 

^ http://www.cfa.harvard.edu/supernova/SNarchive.html 
® However, we highly encourage SN observers to remove telluric 
absorption in their sp ectra as part of their re duction process, using 
standard procedures (JMatheson et al.|2000|l. 

® http://people.lam.fr/blondin.stephane/software/snid/ 


responds to around E-band maximum light, though the 
exact date of maximum is not known. We checked this 
ID by using SNID with the afore-mentioned augmented 
set of SNID templates. With a phase limit of only spec¬ 
tra with phases less than 10 days after maximum, SNID 
finds good matches with both SNe Ic-bl (SNe 2007ce and 
1997ef) and with SNe Ic (1990aa, 2004aw) - though none 
can fit the feature at around 8200 A well (observed wave¬ 


length). Thus, we call SN 2013ez/GRB 130215A a ”SN 
Ic/Ic-bl” in the table, since SNID does not find an unam¬ 
biguous type, and we do not include it when generating 
mean spectra (see below). In addition, since the red- 
shift of this SN/GRB is uncertain and there is only a 
lower limit of 2=0.597 from metal absorptio n lines su¬ 


perimposed on the GRB afterglow spectrum (Cucchiara 


& Fumagalli 2013), the derived velocities based on the 


various absorption lines will be systematically uncertain 
and thus are not included in the velocity analysis below. 

For SNe 2004dn, 2005kl and 2007cl, which have either 
very early or late data where He lines may be too strongly 
visible, we also checked by comparing wi th mean spectra 
of SNe Ib and of SNe Ic for those epochs ( Liu et al.|2015 ) 
and found that they were all clearly not consistent with 
SNe Ib, but rather with SNe Ic. Thus we can call them 
bona-fide SNe I(E3 

In total, we have included 119 spectra of 11 SN-GRBs, 
as well as 120 spectra of 10 SN Ic-bl without GRBs and 
209 spectra of 17 SN Ic, from the l iterature and as par t 


of the Stripped SN data release (Modjaz et al. 
Out of those 448 spectra, 41 were taken later t 


2014). 
ran 90 


days after E-band maximum and thus were taken during 
the nebular phase, where the radiative transfer physics is 
different than during the photospheric phase. Thus, they 
are not analyzed here, but are listed for completeness. 


2.2. Determining the date of maximum light in an 
uniform fashion 

Here, we describe our approach for ensuring that we 
consistently use the date of maximum light in the same 
rest-frame filter as our reference point for all SN e in 
our sample. Following Blondin & Tonry (2007), we 
choose the E-band since most current SIN photometry 
data sets have best-samp led light curves in the E-band 
(e.g.. Bianco et al.||2014 ). For SNe, where the published 
literature had determined their date of maximum light 
in a filter other than the E-band, but published their 
E-band data (i.e., SNe 2004ge, 2011bm), we have deter¬ 
mined the date of maximum in the E-ba nd, using the 
Monte Carlo polynomial fitting method of [Bianco et al.| 
(2014). For high-redshift objects (such as a few of the 
SiN Ic-bl connected with GRBs), we also ensured that 
the filter we used corresponded to rest frame E-band 
filter. For objects where only i?-band or H-band data 
were available (marked with a f in the tabl e, so 8 out of 
34 obj ects), we used the transformations of Bianco et al. 
(2014| (namely [Date of E-band max] = [Date of A-band 
max] - 1.8 days and [Date of E-band max] = [Date of B- 
band max] -|- 2.3 days, respectively, both with a standard 


In |]V[odjaz et al. |2014| we noted that for these SNe we could 
not exclude the possibility that they may have developed or had 
He lines over the phases not covered by our data. Howeve r, with 
the avera ge spectra of SNe Ib and SNe Ic we constructed in|Liu et| 
[aI]pM^ , we can exclude that they are SNe Ib. 
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deviation of 1.3 days, see their Table 10) to convert from 
the R- or il-band date of maximum, respectively, to that 
in the 1^-band. We note that for SN 201 Ibm, the date 
of maximum in the l^-band is different by 11 days from 
that in the i?-band - which is mu ch larger than foun d 


for a large set of Stripped SNe in Bianco et al. (2014). 
However, SN 201 Ibm appears to be an odd SIN Ic, as it 
exhibits large ejecta masses (7—17 Mq), whereas most 


(e.g, Drout et al. 

2011 Cano 

2013 

Lyman et al. 

2014 

bianco et al., m prep). In add 

ition, 

or SiN 20U4te, there 


were multiple datasets in the literature, and we used the 
date o f H-band maximum as measur ed by [Bia nco et ah] 
(2014), which is different from that in Drout et al. (2011|) 
by 1.3 days. 

3. QUANTIFYING SN DIVERSITY BY CONSTRUCTING 
AVERAGE AND STANDARD DEVIATION SPECTRA 

For subsequent analysis, it is important that we remove 
the effect of any host galaxy reddening on the spectra we 
compare, and that we find a robust way to compare spec¬ 
tra taken with a heterogeneous set of telescopes (with 
different binning and wavelength ranges). To that aim 
w e follow the same steps for pre-processing all spectra as 
Liu et al. (2015|)(see details in their Section 3.1), who 


in turn follow Blondin fc Tonry (2007). While the pri- 


mary purpose ot the spectra m 


Blondin &: Tonry (2007) 
library ot SiNlU, a code 


is to be included in the spectra 
for SN classification, here we use them to characterize 
the different SN classes. In summary: we binned the 
spectra to a common logarithmic wavelength scale and 
in that way mapped them onto the common wavelength 
range of 2500—10000 A. We then removed the shape of 
the pseudo-continuum, in order to focus on the charac¬ 
teristic absorption lines, and to not have to worry about 
reddening effects on the continuum (reddi ng is investi¬ 
gated separately in photometry papers, e.g., [Drout et al. 


20111, or about any potential GRB afterglow light (which 


Matheson et al. 

2003 

[Deng et al. 

2005). In 

we checked that the detauit spline 

rom SiN ID 


was reasonable even for the broad spectral features of 
broad-lined SNe Ic. Thus, the rest of the analysis and 
comparisons between SNe Ic and SNe Ic-bl should be ro¬ 
bust as they are based on spectra flattened with the same 
number of spline knots. 

We have already released the co nstruc ted set o f SNIP 
spectral temp lates based on the |M14| data in [Liu &| 
jModjazj (|2014|) , which are publicly accessible via http: 
77wwwTcosmo.nyu.edu/SNYU/ - the newly constructed 
set of SNID spectral templates based on literature data, 
as presented here, are released on our SNYU github 
page (https://github.com/nyusngroup/SESNtemple), 
along with the data products discussed below (in partic¬ 
ular, the mean and standard deviation spectra). 

We note that the maverick PT F12gzk displays very 


high Doppler absorption velocities (Ben-Ami et al.||2012 
and see our Fig. O similar to th ose ot SJN Ic-bl lan d 


also based on radio observations Horesh et al. 2013), 
but does not display broad lines. 'I'hus, we do not in¬ 
clude PTF12gzk when constructing mean spectra (see 
Section]^ for details). 

3.1. What is the “typical SN Ic” or “typical SN Ic-bl” ? 


In a number of crucial situations, one frequently en¬ 
counters the question “What does a typical SN Ic look 
like?” or “Is there a typical SN Ic-bl ?” Such situa¬ 
tions include when intending to determine the type of 
a newly discovered SN, when trying to determine if a 


new transie nt is a novel type of explosion (e.g., Kasliwal 


tion 0 

t SiNe Ic in SNe la searches (e.g., [Poznanski et al. 

2002 

Graur & Maoz 

2013 

Graur et al. 

2015). Histori- 


taken as the best and representative example, as in the 
case of SN 19941 (e.g., Filippenko et al.|1995 Richmond 


et al. 1996), which was a SJN Ic in the nearby iV151, or 
SiN 1998bw/GRB 980425 (Galama et al. 1998) referred 
to as the “typical” SN-GRB, which was the first, and 
closest, instance of the SN-GRB connection. 

However, the natural next step, given the large amount 
of data we have, is to employ a statistical approach, by 
computing the mean and standard deviation spectra of 
a specific type of SNe at a specific phase, to investigate 
the amount of homogeneity, or the lack thereof, in all 
SNe of the same class. We have done so for the spectra 
of SNe Ic and SNe Ic-b l (this work) an d for the spectra 
of SNe Ib and SNe Hb ( Liu et al.|2015 ), the four general 


categories of stripped SiNe . 

We construct the mean spectrum, as well as the stan- 


Blondin et al. 

(20121) and 

Liu et al. 

(2015) 


two specific phases: at maximum light and at 10 aays 
after maximum light (for mean spectra at other phases, 
see Appendix), where all phases are in the rest-frame of 
the objects. For all cases, we included spectra that are 
within ±2 days of those dates. On Fig as well as on 
subsequent figures, we also note the number of individual 
spectra that were used to construct the mean spectrum 
for each epoch range. Sometimes, there is more than one 
spectrum per SN included (especially for mean spectra 
at maximum light) as we included multiple spectra of 
the same SN from various groups with different wave¬ 
length ranges, as well as multiple spectra of the same 
SN obtained over the four-day bin. In order to check 
that the mean spectra are not dominated by a few SN 
with large numbers of spectra, we also constructed mean 
spectra, where we include only one spectrum of each in¬ 
dividual SN. There is no statistical difference between 
mean spectra based on one spectrum per SN vs. those 
based on multiple spectra per SN. We note that the mean 
spectrum of broad-lined SNe Ic include more spectra of 
SN-GRBs than those of SN Ic-bl without observed GRBs 


(see Section 3.3), such that it is weighted more towards 


SN-GRBs spectra. This is due to the apparent higher in¬ 
terest in observing and following up SN-GRBs than SN 
Ic-bl without observed GRBs. 

In addition, we test whether the so-called “prototypi¬ 
cal” SN Ic 19941 and SN Ic-bl SN 1998bw are truly rep¬ 
resentative of the respective SN sub-types by comparing 
them to the mean spectra for that class. Fig. shows 
that SN 19941 is more than one standard deviation away 
from the mean SN Ic spectrum, even reaching the maxi¬ 
mum deviation at many wavelengths, meaning its spectra 
are the most deviant of the SN spectra in the sample - 
at both maximum light and 10 days after max. Thus, we 
conclude that SN 19941 cannot be called a typical SN Ic 
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Fig. 1.— Mean spectra of SNe Ic (in black, left) and of all SNe Ic-bl (including both with and without GRBs, right) with corresponding 
standard deviation (in yellow) and maximum deviation (in orange), after we removed the pseudo continuum (see text for more details), at 
peak brightness di 2 days (top), and at 10 days after peak ± 2 days (bottom) in the rest frame. SNe 2013ez and PTF12gzk are not included 
here (see text). Shown are also the number of individual spectra th at were used to construct the mean spectrum for each epoch range. For 
comparison, overlai d is the SNID-ified spectrum of SN 19941 (left, [Filippenko et al.||1995| |Modjaz et aL||2014|l and of SN 1998bw (right, 
|Galama et al.[|1998|l, which have been considered as the standard SN ic and bis] ic-bl, respectively, me spectra of SN 19941 are more than 
one standard deviation away from the mean spectra of SNe Ic, often reaching maximum deviation, at many wavelengths. On the other 
hand, SN 1998bw appears to be a typical SN Ic-bl - its spectrum is very close to the mean spectrum and well within one standard deviation, 
specially at 10 days after max. We note that there were twice as many SN-GRBs than SN Ic-bl without observed GRBs that went into 
constructing the SN Ic-bl mean spectrum, and thus, that our mean SN Ic-bl spectrum is weighted more towards SN-GRBs spectra. Note 


that the apparent line at 7400—7800A and its large standard deviation in the mean spectra is prob ably d ue to a combination of O I A7774 
and telluric (i.e., atmospheric) absorption (A-band at 7620 A); the latter was removed only for the M14|and Berkeley spectra. 


spectrally, in contrast to common belief. The comparison 
of its photometric properties reveals a similar and even 
more pronounced conclusion, na mely that it is not repre- 


From the massive stellar evolutionary point of view, it 
is very hard to produce a SN-GRB progenitor which is 
both H- and He-free at the low metallicites at which most 


sentative of light c urves of SNe Ic (Richardson et al. 2006 
Drout et al. 2011 F. Bianco et al, m prep.), since it fadec 


much faster than the average SN Ic light curve. On the 
other hand, SN 1998bw seems to be a “typical” SN Ic-bl 
(Fig.0 - at both phases, SN 1998bw is consistent with 
the mean spectrum well within one standard deviation 
at most wavelengths for its spectrum at t = +0 days, 


SN-GRBs are observed (Modjaz et al. 2008 Levesque 


the n 

lodels tor GRB progenitors are He-stars (Woosley 

et al. 

2002 

Yoon & Langer|20051|Woosley & Heger|20061 

Yoon 


, with He masses between I.O — 15 Mq m the 


Age-Main-S equence masses of 15 — 50 Mq (Woosley 
et al.|[2002). While hiding He due to insufficient mixing 


days. While we only show mean spectra at two specific 

has been considered iHachinger et al. 12012 

Dessart et al. 

epochs, our conclusions remain unchanged when making 

2012 

Biro & iVlorozova|2014 Liu et al. 2015 

, the amount 

this comparison at other phases. Thus, using SN 1998bw 

that can be hidden appears to be very small (<U.I iVirr^ 

as a spectral template for other broad-lined SN Ic and 

of He; 

Baron et al. 

1996[ 

Hachinger et al. 

2012 

), much 


warranted, but note that our SiN Ic-bl mean spectrum is 
more weighted towards spectra of SN-GRBs. 

3.2. The Helium Problem for SN-GRBs - is there 
Helium hidden in SN Ic-bl? 


smaller than the expected mass of the He progenitor star. 
Another possibility may be that the He lines are present 
but not detectable in SNe Ic-bl spectra due to smearing 
by the high ejecta velocities in SNe Ic-bl. The latter is 
what we want to address here by simulating the smear¬ 
ing due to high velocity on real SNe Ib spectra, i.e. those 
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with detected He lines - effectively simulating a “broad 
lined SN Ib”. We then want to explore if those broad¬ 
ened He lines could have been detected in spectra of SNe 
Ic-bl and SN-GRBs. 
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based on the results of our template fitting method dis¬ 
cussed below (see Fig.[^. We note that in order to match 
the spectral features of the mean SN Ic-bl spectrum, in 
addition, we blue-shifted both the convolved mean SN 
Ib and the mean SN Ic spectra by 3,000 km s“^. A 
convolution with a Gaussian profile (or inverse Gaussian 
smoothing) simulates well the smearing of th e spectra 
due to incr ease of kinetic energy per unit mass (|Iwamoto 
et al.|1998| K. Nomoto, private communication), i'he re¬ 
sults are shown in Fig. while convolving a SN Ic mean 
spectrum with a Gaussian profile yields a spectrum that 
is very similar to a SN Ic-bl mean spectrum (except over 
a small wavelength range of 4000—5000 A, most likely 
due to Fe blending, see below), a similarly convolved SN 
Ib mean spectrum is different from the observed SN Ic-bl 
mean spectrum, by more than 3-4 standard deviations at 
the positions where the He lines are expected, especially 
for spectra at 10 days after maximum. The smeared- 
out He lines in the convolved SN Ib spectrum are even 
outside the maximum deviation spectra of the observed 
SN Ic-bl mean spectrum, especially for spectra days after 
maximum light. The match between the SN Ic-bl and our 
smeared SN Ic spectrum is not unique; a small range of 
convolution kernels and blueshifts provide similarly good 
fits. In all such cases, using the same fit parameters for 
simulating a “broad-lined SN Ib”, the broadened He-lines 
should have been detected in spectra of observed SNe Ic- 
bl, including those of the SN-GRBs, but they are not. 
Indeed, the choice for the convolution parameter shown 
in Fig.l^is more than conservative, since it is at the high 
end of the range of convolution parameters when doing 
a formal template fitting procedure for the Fe II A5169 
line, meaning it is amongst the most extreme, i.e. worst, 
cases of smearing-out for SNe Ic-bl and SN-GRBs (see 
Fig.[§). 

Our conclusions remain unchanged when performing 
the same comparisons for mean spectra at different 
phases, as long as they are after maximum light, when 


Rest Wavelength [A] 

the He signature is strong enough in SNe Ib for a robust 

detection ( 

Liu et al. 

2015 

). In the same spirit, we re- 

Fig. 2.— Mean SN Ic-bl spectra for two different phases (in 

peated this exercise 

Dy simply adding He lines into the 


black), as well as an average SN Ic and an average SN Ib spec¬ 
trum - the latter two convolved with a Gaussian and blue-shifted. 
The average SN Ic-bl spectrum can be reasonably reproduced by a 
normal SN Ic spectrum convolved with a Gaussian with (7=6,000 
km s“^Gaussian (i.e., FWHM=14,000 km s“^), and blueshifted by 
3,000 km s“^. However, a similarly convolved SN Ib spectrum de¬ 
viates from the average SN Ic-bl spectrum at the positions where 
the He-lines are expected, and is well outside the standard and 
maximum deviations (the latter not shown) - indicating that even 
a broad-lined SN Ib, i.e., a smeared out SN Ib spectrum, would 
have had detectable signs of He lines, which are not observed in 
our observed SN Ic-bl mean spectra. Annoted are the lines of 
He I A5876, He I A6678 (which may be blended with Ha, a s SNe 
Ib may show some Hq; |Parrent et al.||2016l |Liu et al.||2015^ , and 
He I A7065 in the spectrum ot the simulated "biN ib-bl '. 


For t his p urpose, we take the SN Ic mean spectra from 


SN Ic mean spectra. Again, we find that smeared-out He 
lines do not match the appropriate region in the SN Ic-bl 
spectrum, even when using very high smearing velocities 
(up to 25,000 km s“^ with respect to the SN Ic spec¬ 
trum). Thus, we conclude that even a smeared-out SN Ib 
spectrum (i.e., a simulated SN “Ib-bl”), would have had 
detectable He lines, which are not observed in spectra of 
SNe Ic-bl and SN-GRBs, nor is it possible to smear out 
the He lines sufficiently make them non-detectable. The 
lack of He lines in SN Ic-bl-GRB spectra is i n line with 


findings by e .g., Taubenberger et al. (2006) and Ghornock 


et al. (2010), who obtained INIR spectra of SNe Ic-bl and 


Section 


3.1 


as well as the S N Ib mean spect rap] (con- 
structeci in'the same way, see Liu et al. 2015) and con¬ 
volve both with a Gaussian of width cr=6,000 km s“^(i.e., 
FWHM=14,000 km s“^), choosing a conservative value, 

The SN Ib mean spectrum for maximum light is based on 38 
spectra of 13 SNe Ib, while the SN Ib mean spectrum for 10 days 
after maximum is based on 37 spectra of 18 SNe Ib. 


did riot observe the stro ng l^m line of He. Recently, 
n (|2015|) claimed a potential detection 
it He 


Milisavljevic et al 
ot the 1 fiia line of He in SN Ic-bl 2012ap via SYNOW. 
However, since they do not detect the otherwise clean 
and strong 2 /im line of Helium in the same spectra at 
phases where one expects it, their evidence for significant 
He detection is weak and the l/im line feature may be 


SNe ( 

Millard et al.|1999 

Gerardy et al. 12004 

Sauer et al. 

2006 

h'riesen et al. 112011 

I. 
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Fig. 3. — As in Fig.ln but for average spectra of SN Ic-bl without observed GRBs (left) and SN Ic-bl with GRBs (right), for two different 
phases. SN-GRB haveoroader lines than SN Ic-bl without GRBs (see text for more details). 


Thus, we conclude that the absence of clear He lines 
in optical spectra of SNe Ic-bl is not due to the smearing 
out of the Helium lines by the high velocities present, as 
otherwise they would have been still detectable, assum¬ 
ing that the same process that operates on all the other 
lines affects the Helium lines in the same way. We discuss 
the implications these findings have for the understand¬ 
ing of SN-GRB progenitor systems in Section Never¬ 
theless, we note that there is a large discrepancybetween 
the simulated SN Ic-bl spectrum (i.e., after broadening 
the SN Ic mean spectrum) and the real SN Ic-bl mean 
spectrum in the blue region (3500-4400 A, Fig.l^, where 
the forest of Fe lines lie, which may have implications 
for understanding line formation in highly expanding at- 


spectra of SNe Ic-bl with GRBs and plot them for two 
epochs (around maximum light and 10 days after) in 
Fig. Inside- by-side. We include SN 2009bb amongst the 
SN-GRBs, since it was suggested to be a relativistic SN 
based on its large radio-emission, similar to nearby SN- 
GRBs, even if no gamma-rays were detected (Soderberg 

:|[20l0t . 

Figure 


et al. 


sis models such as CFMGEN 

Dessart & Hillier 

2007), 

PHOENIX iBaron et al. 

1995 

, and SEDONA 

Kasen 


been done before. 

3.3. Differences in Spectra between SNe Ic-bl with and 
without GRBs 

Two of the outstanding questions in the field are what 
distinguishes SNe Ic-bl with observed GRBs from those 
without observed GRBs, and what produces SN Ic-bl, 
including those without observed GRBs. While host 
galaxy analysis have shown that t here may be differ¬ 
ences in metallicity and SFR density d Modjaz et al 12008 


Sanders et al. 2012a Kelly et al. 2014p, here we study the 

explosion properties, and in particular the early-time, 
photospheric spectra. To that aim, we construct mean 
spectra of SNe Ic-bl without observed GRBs and mean 


in Figure we can see that the features in SNe Ic-bl 
with GRBs are somewhat broader than those in spectra 
of SNe Ic-bl without GRBs, especially at 10 days after 
R-band maximum. For example, for the SNe Ic-bl with¬ 
out GRBs at tvmax = +10 days, the wavelength range 
between 4000A and 5000 A appears to have a number of 
’’features” — 2-3 troughs —, while the same wavelength 
range for the SN-GRBs shows only very smooth undula¬ 
tions and at most one broad feature. We note that our 
conclusions remain unchanged even if we restrict our¬ 
selves to include only one spectrum per included SN, in 
order to avoid that many spectra of the same SN domi¬ 
nate the mean spectra. As shown below in Section |4.4[ 
this trend is also seen quantitatively in the amount by 
which the Fe II A5169 line in the mean spectrum of SNe 
Ic has to be broadened in order to match the SN Ic- 
bl spectra (i.e. the convolution velocity of the template 
fitting method, see also Appendix Section 0. We will 
address this question using a complementary method, 
namely measuring the Doppler-shifted absorption veloc¬ 
ity of Fe II A5169 in Section where we find the same 
trend, but for Doppler-shiftea absorption velocities: SN 
Ic-bl/GRB have higher velocities, on average, than SNe 
Ic-bl without observed GRBs. We will discuss the impli¬ 
cations of these observations for the SN-GRB explosion 
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mechanisms and production mechanisms of GRBs in Sec¬ 
tion 0 

4. SPECTRAL LINE VELOCITY EVOLUTION 

Now we shall shift our attention to measuring veloci¬ 
ties, since they provide clues about the dynamics of the 
explosion, and since an estimate of the photospheric ve¬ 
locity is needed for computing important quantities, such 
as ejecta masses, from light curves and spectra (e.g., 
Drout et al.|201l]|Caiio|2013||Lyn:ian et al.|2014[ l^ddial 

et al.|2015p , using analytical models (e.g., |Arnett| ( |1982p . 

Here, it is crucial to distinguish between two main kinds 
of velocities: the absorption velocities, called Vabs, which 
are measured from the (Doppler) blue-shifted trough of 
specific lines, and the width velocities, by which the lines 
are broadened (which we call VconvoiuUon based on the 
method we use, see Appendix®. We measure both sep¬ 
arately — show in Section |4.d| that different lines and 
different methods give different absorption velocitie s — 
and present their self-consistent analysis in Section 
for Fe II A5169 absorption velocities, and in Section 


4.4 


4.6 


for the width velocities respectively. As we will see below, 
while for most cases these two kinds of velocities corre¬ 
late with each other and are very similar, as expected 
from radiative transfer models in t he kinds of expandin g 
atmospheres considered here (e.g., |Dessart et al.||20Il|), 
there are outliers: cases where a SIN has high absorption 
line velocities, but small line widths (e.g., PTF 12gzk) 
and vice versa, a SN that has large line widths but low 
absorption velocities (e.g., PTF lOvgv). Thus, until a 
unifying picture of spectral line formation in SN Ic and 
SN Ic-bl is developed and tested (e.g., including aspheric- 
ities in the ejecta geometry), we suggest that authors, 
when reporting velocities, are clear about how the num¬ 
ber was derived (whether from modeling or observations, 
and which line was used) and whether it is the absorption 
velocity or the width of the line that they are reporting. 

4.1. Absorption line velocities and their measurement 

errors 

We determine the b lueshi fts in the troughs of the ab¬ 
sorption lines from the |M14| Stripped SNe sample, as well 
as those from the literature (see Table 1). For the SNe Ic 
spectra, we use the robust technique of measuring SN line 
profiles and of e stimating error bars as presented by us in 
Liu et al. (2015). Here we summarize the salient features: 
we constructed noise-smoothed spectra and measured the 
location of ma ximum absorpt i on usi ng the procedure rec¬ 
ommended by Blondin et al. (20121 of fitting a parabola 
to a small region around the minimum absorption. It 
is especially powerful since it does not assume specific, 
Gaussian, for the full line profile and therefore allows for 
asymmetric shapes for the full line profile. We refer to 
our velocity measurements as absorption line velocities 
(in short ”vabs ”) in the rest of this work. We verified that 
our velocity measurements are robust: they yielded the 
same values as on the continuum-removed spectra, which 
were produced during the const ruction of mean spectra. 
As we showed in Liu et al. (20151, we constructed spectral 
uncertainty arrays to be used tor the subsequent estima¬ 
tion of the velocity measurements errors. We constructed 
spectral uncertainty arrays from the actual spectra them¬ 
selves by separating noise from signal in Fourier space, 
taking advantage of the fact that the typical width of 


a supernova spectral feature is large, since it is Doppler 
broadened by the ~ 5,000 - 20,000 km s~^ expansion ve¬ 
locity of the SN ejecta. For computing the errors on the 
line identification velocity measurements, we conducted 
Monte Garlo simulations where we generated 3,000 mock 
spectra, using the associated spectral uncertainty arrays, 
and reran our automated velocity measure ment code on 
them (for more details see|Liu et al.||2015). 

We did this analysis for Fe H A5169 , since Fe II A5169 
has be en suggested to tra ce well t he photospheric ve- 


locity (Branch et al. 2002 but see Dessart et al. 2015 
for caveats). On the other hand, as we discuss bkow, 
while Si II A6355 is commonly identified with the most 
isolated feature at ~6000—6100 A in broad-lined SNe Ic 


(e.g.. 

Patat et al. 

2001 

Chornock et al. 2010 

Lyman 

et al.| 

2014|), the identification ot Si II A6355 in be 

)th SJNe 


et al.|20f6 1 . We note that our line identifications are not 


definitive since we have not conducted spectral synthesis 
calculations for all spectra. For clarity, in Appendix 
we describe how we self-consistently identify Fe II A5169 
for all subtypes of SN Ic. 

However, this procedure (which we call the ’’line iden¬ 
tification method” hereafter) does not work properly for 
spectra of SNe Ic-bl, since those spectra are severely 
blended - in particular, the multiple lines of Fe II 
(Fe II A4924 , Fe II A50I8 , and Fe II A5169 ), which are 
visible in SNe Ic, are blended into one broad feature in 
SNe Ic-bl (see Fig. [^. As we detail in Appendix [A] here 
we develop a novel method to measure line velocities and 
widths for SNe Ic-bl using a template-fitting method. In 
summary, we use the mean spectrum of SNe Ic that we 
constructed in Section^ (in particular the region around 
Fe II A5I69), blue-shirt it and convolve it w ith a Gaus¬ 
sian (similar to what we did in Section 3.21 to find the 
convolution width and blueshift that best matches the 
input SN Ic-bl feature. We implemented this process by 
using the Monte Garlo Markov Ghain (MCMC) sampler 
’’emcee” (Foreman-Mackey et al. 20I3|) to find the best fit 
and to obtain error bars, 'i'his procedure usually yielded 
very good fits as indicated by the ^ unity values for the 
min values (e.g., see figure 10). We tested that both 
methods (the identification method and the template fit¬ 
ting MGMG method) yielded the same values of Vabs for 
SNe Ic for which both methods can be applied. Thus, 
we think that there is no systematic shift or bias when 
comparing the Vabs values of SNe Ic (obtained via the line 
identification method) to those of SNe Ic-bl (obtained via 
the convolution MCMC method). We compare our de- 
rived velocities to those from the literature in Section |4^ 
For the low S/N spectra of SN 2010ma/GRB I0I2I9B, no 
satisfactory fits could be achieved, so we do not include 
any of the fitting results for SN 2010ma/GRB 10I2I9B. 

We emphasize that even if the feature we identify as 
Fe II A5I69 is not due to Fe II A5I69, relative com¬ 
parisons between the different SN types will still hold 
and are valuable, since we are self-consistently using the 
same feature for velocity measurements for all SN Ic sub- 
types. As to Si II A6355, there are some suggestions 
that the commonly identified feature may not be due to 
Si II A6355, but rather a comb ination of elements, in¬ 
cluding Ha ( Parrent et al.|20T6 ). Furthermore, there are 
sometimes two features (one at 6100 A and one at 6300 
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TABLE 3 

Measured absorption 
VELOCITIES 


Phase^ 

(days) 

VabsfFe II A5169) 
(km s“^) 

-15.0 

snl983V 

-18000 ± 970 

-12.0 

-19000 ± 1100 

-10.0 

-19000 ± 880 

-2.0 

-12000 ± 940 

5.9 

-10000 ± 460 

6.9 

-10000 ± 400 

7.9 

-10000 ± 430 

8.9 

-9900 ± 540 

9.9 

-9500 ± 530 

10.9 

-9300 ± 510 

11.9 

-9100 ± 550 

31.8 

-9000 ± 680 


Note. — This table is avail¬ 
able in its entirety in a machine- 
readable form in the online 
journal. A portion is shown 
here for guidance regarding its 
form and content. 

^ Rest-frame age of spectrum in 
days relative to E-band maxi¬ 
mum. See text for details. 


A) in that range, and without spectral synthesis calcula¬ 
tions for all spectra it is hard to identify which one is due 
to Si II A6355. While we have measured the absorption 
wavelength of the strong feature at ^6300 A for SNe Ic, 
and for SN 19941, we show absorption velocities based on 
the identification of that feature with Si II A6355 only for 
comparison purposes (Fig.|^ and do not report them for 
our full sample. In addition, for SNe Ic-bl we find from 
our template-fitting that the SN Ic template (after blue 
shifting and broadening) is not a good representation of 
the SN Ic-bl spectra in the wavelength range of ~5800 
A— 6300 A and we do not obtain good fits (Appendix [ a]). 

Using the two different methods above, we list the mea¬ 
sured velocities and their error bars for the Fe II A5169 
line as a function of spectral phase for all SNe in our 
sample in Table This table is available in its entirety 
in a machine-readable form from the online journal. The 
listed and plotted error bars represent the 68% uncer¬ 
tainty range. 


4.2. Comparison with literature values of Absorption 
velocities for SN-GRBs and Stripped SNe 

Two recent papers have measured velocities in large 
samples of Stripped SN spectra that include SNe that 
are analyzed in this work. Since they did not compute 
velocity uncertainties, the following comparison will in- 
cl ude only our velocity error values. 


Schulze et al. (2014) include the same seven SN-GRBs 
as we, mostly with the same published data. However, 
they attempt to measure Fe II A5169 with a method 
that is different from our main method of template fit¬ 
ting, but somewhat similar to our line identification 
method: they attribute the blended feature at ~ 5000 
A to Fe II A5169 alone. Of the seven SN-GRB, where 


we and Schulze et al. (2014) use the same spectra pub¬ 
lished by them and others, tneir values for five SN-GRBs 


are consistent with ours within our error barJ^ For 
SN 2013cq/GRB 130427A, where both velocities are for¬ 
mally consistent within the error bars, their value is 
larger by ~ 6000 kms“^. Also, for SNe 1998bw and 
20031w, their velocity values are all systematically larger 
than ours (by 12,000 km s~^ on average, which corre¬ 
sponds to more than 3-(T errors). This discrepancy is 
most likely due to the fact that their method assumes 
that the whole feature is only due to Fe II A5169, and 
they do not include any contributions from the other 
two Fe II lines (namely Fe II AA 4924, 5018), which are 
blended with Fe II A5169 in many of the spectra of SNe 
Ic-bl (see Fig. |^. In addition, without attempting to 
give an interpretation, we note that for some SN Ic-bl 
spectra the relative strength of the two parts of the Fe II 
” W” feature are the opposite compared to the SN Ic tem¬ 
plate spectrum. This is the case for SN 1998bw and SN 
2002ap at some, but not at all, phases. However, this 
does not affect our velocity measurements. The tem¬ 
plate fitting is most sensitive to the width of the entire 
complex, and in all instances where the relative feature 
strength is opposite to that in the template, we find that 
the edges of the feature coincide with the edges of the 
stretched and shifted template, and we conclude that 
the broadening and location of the feature are properly 
measured. For SN 1998bw, in particular, we think that 
our method is supe rior to the identification method (e.g.. 


Schulze et al.|2014 ), as the latter identifies the wrong fea- 
ture with Fe II A5169, namely the one that is actually due 
to the blend of Fe II A4924 and Fe II A5018. We think 
that this mis-identification is due to the fact that, as 
mentioned above, for SN 1998bw/GRB 980425 the fea¬ 
ture due to Fe II AA 4924, 5018 is stronger than the one 
due to Fe II A5169, which is the ot her way around com- 
pared to most other SNe, and that |Schulze et al. (2014) 
presumably used the dominant feature in that complex 
to identify as Fe II A5169. We did try to perform a fit in 
which we forced the dominant feature at ~4700A to be 
due to Fe II A5169 but those fits were ruled out. 

We have 14 SNe of Types Ic and Ic-b l in common with 
the SN sample of Lyman et al. (2014) who use mostly 
the same data sets as we do, but measure Fe H velocities 
for some SNe and Si II A6355 velocities for others (from 
spectra taken closest to maximum light) in a highly inho- 
mogenous way. While we fit a cubic spline to the distinct 
absorption feature of Fe II A5169 in SN Ic spectra and 
fit a SN Ic template (consisting of the SN Ic mean spec¬ 
trum) to the blended absorption features of triple Fe H in 
SN Ic-bl spectra, they fit a simple Gaussian profile to ab¬ 
sorption features of Fe H or Si H. However, for some SNe, 
they take values from the literature (either from spectral 
synthesis or from interpolating the velocity values from 
spectra taken at epochs far from maximum light) and 
for others, even use mean values within SN subtypes in 
cases they could not identify Fe or S|^ For seven of the 
SNe, the velocity values are consistent with ours within 


For SN 2003d h/GRB030329. we no te that its phases are not 
correct (Fig. 5 of [Schulze et al.]|2014[ |, as the phases were de- 
redshifted twice (G. Eelouclas, private communication). Adopting 
the correct phases, their values are consistent with our measured 
values within the uncertainties. 

We note that the total numbe r of SNe for all types we have in 
common with [Lyman et ah ]( |2014| l is 26 and the conclusions using 
the full data set remain tne same. 
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our error bars (ie. within 1 a, while for four SNe we 
are consistent at 1—2 a and three SNe we are co nsistent 
within ~2 a . Given that the measurements of |Lyman| 


et al. (2014) are mostly on the same spectra as ours, we 
can sirnply assume that their velocity measurements have 
the same uncertainty; in that case their agreement with 
our values is even stronger, such that all their measure¬ 
ments agree with our values within each measurement’s 
uncertainty. 

4.3. Different methods yield different velocities 



Fig. 4.— Velocity evolution of different lines and using differ¬ 
ent methods for the same SN, SN 19941, using the same data, 
in almost all cases. The absorption (i.e., Doppler) velocities val¬ 
ues for Si II A6355 and Fe II A5169 were measure d by us, while 
the SY NOW photospheric velocities are taken from [Millard et al.| 
l|1999|l and the velocities as computed from the M onte Carlo SJN 
spectrum synthesis co de ofiMazzal i et al .l ||20 08|| (and referenc e 
therein) are taken from|Sauer et ar] l|2UUb|rancl|r'ian et all] (|2006||. 
Note that different methods give ditterent values, with bVlMUW- 
based ahd Fe II A5169 velocities systematically higher than Monte- 
Carlo-based and Si II A6355 velocities. We note that there is some 
controversy whether the feature at 6100—6200 A is purely due to 
Si II A6355 ( [Parrent et al.|2016[ |. 


Since it is not widely known that different absorption 
lines yield different Doppler veloci ties because of optica l 
depth and ionization effects (e.g., Mazzali et al. 2013), 
here we use SN 19941 to show the range and types of dif- 
ferences, since this SN was extensively analyzed by dif¬ 
ferent group s, including theoretical m odeling by various 
groups (e^., Baron et al.||1996 1999 and see below). 

In Fig. we plot two velocities computed from mod¬ 
eling in the literature, and absorption Doppler velocities 
measured by us for two different lines in SN 19941, all 
based on the same or similar data: photospheric expan¬ 
sion velocity as comp uted from the SN Monte Carlo spec- 


trum synthesis code (Mazzali 2000 Mazzali et al. 2008 
and refe rences therein) as applied to SIN 19941 (as pre¬ 
sented iniSauer et al.|2006[| Pian et al.|2006|), the SYNOW 
photospheric spectra frorn Millard et al. 1 1999), and the 
absorption velocities for Fe II A51t)9 and Si II A6355, as 
measured by us. In terms of the various spectral synthe¬ 
sis codes, different codes have different degrees of simpli¬ 
fications and assumptions, which is beyond the scope of 
this paper to describe in detail. Broadly speaking, cer¬ 
tain codes such as the simple SYNOW code assumes a 


sharp photosphere, though SYNOW includes accommo¬ 
dating ’’detached” layers for certain elements (e.g. H or 
He) that may move at higher velocities. 

The data we used for m easuring the two absorp tion 


lines, namely the Berkeley (Filippenko et al. 1995) and 
|M14| spectra, constitute all or almost all of the same 
d ata as used for mo deling in Millard et al. (1999) and 


Sauer et al. (2006). The highest velocities come from 
the SYNOW models, as well as from the absorptions 
velociti es of Fe II A5169 , c onsistent with the sugges¬ 
tion of |Branch et al.| (|2002|) that, based on SYNOW, 
the Fe 11 A5169 line is a good tracer of the photospheric 
velocity of the ejecta. The MC photospheric velocities 
are lower than the aforementioned two - and impor¬ 
tantly, they have a different velocity evolution shape, as 
they do not plateau at late times like the SYNOW and 
Fe 11 A5169 velocities, but continue to decrease rapidly. 
Thus, at later times (A tymax > +26 days) the difference 
is up to 7000 km s“^between MC velocities and both 
SYNOW and Fe II A5169 velocities. The Si II A6355 
velocity behavior is similar to that of the MC velocities, 
but has a slight offset to lower values (by ~2000 km s“^). 
Thus, velocities based on SYNOW code and Fe II A5169 
measurements have a similar evolution on the one hand, 
while on the other hand, those based on the MC code and 
Si II A6355 m easurements have ano ther evolution. This 
explains why Tanaka et al. ([2009), for example, claim 
that for SN 2008D, the Si 11 A6355 line (and not the 
Fe II A5169 line) is tracing reasonably closely the posi¬ 
tion of the photosphere, since they compute photospheric 
velocities using the MC code, not SYNOW. 

The offset between SN 1994rs MC photospheric veloc¬ 
ities and measured Si II A6355 velocities (such that the 
MC photospheric velocities are higher by ^2000 km s“^ 
than the Si II A6 355 velocities) was already noted in 


Sauer et al. (2006) (their Fig. 5). They suggested that 
the velocity offset could be due to two reasons: the line 
at ~6100 A (assumed be solely due to Si II A6355) may 
have significant contribution from Ne I 6402, such that 
the measurements do not trace Si II A6355 well; or there 
is a real discrepancy between measured Si II A6355 ve¬ 
locities and the computed photospheric velocities. 

While we performed this analysis explicitly only for 
SN 19941, which indeed may not be a typical SN Ic (see 
above), the same findings hold when comparing the pho¬ 
tosph e ric vel ocities computed by the MC code of Mazzalij 
et al. ( |20^) (for a raiige of SN-GRBs, as recently pre- 


sented in Walker et al.||2014 and references ther ein) and 


the measure d Fe velocities tor the same SN (e.g^ Schulze 


et al. 2014 and as presented here in Section WTr the 


MC-based velocities are systematically lower than the 
measured Fe absorption velocities. 

While we are not trying to resolve the issue of which 
line to use, nor advocate for using a specific method, we 
do recommend that any meaningful comparison of veloc¬ 
ities between different SNe and different kinds of SNe Ic 
should be done using the same method and the same 
line for velocity determination. So far, modeling has 
been done mostly for the SN-GRBs, while Si 11 A6355 or 
Fe II A5169 measurements have been reported for SN Ic 
in the literature, and only Si II A6355 measurements for 
SN-GRBs since all other lines appear to be too blended. 
Especially Si II A6355 vs. Fe II A5169 velocities should 
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not be used interchangeably, given their large offset from 
each other (up to 7,000 km s~^ at 2—3 weeks after max¬ 
imum) and neither should the ” photospheric” velocities 
derived by SYNOW vs. the MC code, given the system¬ 
atic offset amongst the two. 

4.4. Absorption velocities of SNe Ic, SNe Ic-bl with and 
without GRBs 

We concentrate on measuring Fe II A5169 velocities in 
the same consistent way for all members of the SN Ic fam¬ 
ily: regular SNe Ic, broad-lined SNe Ic without GRBs, 
and SN-GRBs. The goal is to quantify in a consistent 
way by how much the absorption velocities of SNe Ic-bl 
are higher than those of SNe Ic, and if there are system¬ 
atic differences for SNe Ic-bl with and without GRBs. 
Spectral synthesis studies suggest that Fe lines are good 
tracers o f the photospheric velocity, since they do not 
saturate ( Branch et al.||2002 ). 

Figure l^shows the h'e II A5169 absorption velocities of 
11 SN-GRBs, 10 SNe Ic-bl and 17 SNe Ic, color-coded 
by SN subtype. This is the most detailed and consistent 
analysis on the largest data set to date. The left panel 
displays the measurements for the individual SNe, while 
the right panel displays the rolling weighted mean veloc¬ 
ity evolution for each SN type in order to dis play their 
bulk e volution. We follow our procedure from lLm et al. 
(2015) to construct rolling weighted means, in order to 


mitigate the impact of the choice of the beginning and 
end points of the bin by replacing the bin with a rolling 
window. To summarize, we took the weighted average 
within a certain window size and moved, or ’’rolled”, in 
1-days steps. We chose the rolling window size for early 
times (before 30 days after max) to be 5 days, as the ve¬ 
locity evolution is fast at early times, and for phases af¬ 
ter that, we chose a 10-day-window. In order to mitigate 
the impact of a single SN that has many data points over 
the window range, we took the weighted average of the 
data points for the same SN (within the window range) 
first. We then computed the weighted mean absorption 
velocity of all SNe of the same type within that rolling 
window if there were more than three SNe with data in 
that window (for any standard deviation to be meaning¬ 
ful). For the center phases of the rolling window, we took 
the mean of the phases at which the velocity measure¬ 
ments were made that contributed to the weighted mean 
velocity. Obviously, the rolling weighted mean velocities 
are not independent of each other, as adjacent windows 
may include many of the same SN data . 

The following was not encountered in Liu et al. (2015|) 
but is addressed here: for the asymmetric errors of 
Doppler absorption velocities for the SNe Ic-bl (derived 
from the template fitting method via MCMC which out¬ 
puts the full posterior prob ability d e nsity function), we 
followed the prescription of Barlow (2003| to construct 
weighted averages from asymmetric uncertainties (their 
model 1, though there was no significant difference be¬ 
tween using their model 1 vs. their model 2). The er¬ 
ror bars for the weighted mean values are the unbiased 
weighted standard deviation of the data (not that of the 
weighted mean), since we would like the error bar to re¬ 
flect the variation of the velocity evolution in the SNe, 
not just the quality of their velocity measurements (i.e., 
the uncertainty in the velocity measurements). In Ta¬ 
ble 1^ we list the weighted mean velocities for both lines 


at R-band maximum, in a window of five days. 

Now, for the first time, we have a large enough 
sample analyzed in a systematic manner to quantify the 
similarity and behavior of SNe Ic-bl (especially those 
accompanied by GRBs) compared to those of SNe Ic 
- something that has been attempted before but only 
by u sing SN 19941 as a representative of all SNe Ic 


(e.g.,|Pian et al.||2006[ |Valenti et al.||2008| |Taubenberger| 

et al.puQ i ^ . A number of important trends can be seen 

in Fig. [Sp where in the following the quoted weighted 
mean vdocities are computed for R-band maximum in a 
five-day window, and where the error bar is the standard 
deviation in the data: 


• When grouping together all SNe Ic-bl (both with and 
without GRBs) in one category, the weighted mean ve¬ 
locity is < vsNic-bi >= —22, 200 ± 9,400 km s“^ . This 
is more than a factor of two above the value for SNe Ic 
(weighted < v > snic = — 9j 700± I,600j^ We note that 
the absorption velocities of SN 2010bn7GRB 100316D 
are the highest overall for any SN-GRB, higher by ^ 
150% compared to any other SN Ic-bl with a GRB at 
R-band ma ximum. This was already noted by |Schulze| 


et al. (2014), however, their measured values of the other 
SiN-GRBs velocities were also high, since they were usin g 
the more simplistic identification method (Section |4.2[ ), 
which tends to overestimate absorption velocities tor 
when Fe II A5169 is highly blended, as in SNe Ic-bl. In 
case SN 2010bh/GRB 100316D is different from the rest 
of the SNe Ic-bl, we exclude it when we list the weighted 
mean velocity of all SNe Ic-bl in Table We find that 
the weighted mean velocity for all SNe Ic^l is still almost 
a factor of two higher than the weighted mean velocity 
for normal SNe Ic. 

• We find that SNe Ic-bl with GRBs have, on average, 
consistently higher Fe II A5169 absorption velocities 
than SNe Ic-bl without GRBs (by ~3,000 km s~^ for 
the date of R-band maximum, and by ^6,000 km s~^as 
an average over all phase, for both cases when excluding 
SN 2010bh/GRB 100316D, see Tabled). SNe Ic have the 
lowest absorption velocities. Nevertheless, as the plots 
indicate, there is an overlap in absorption velocities be¬ 
tween the different types of SNe, though the systematic 
trend of highest mean velocities for SN-GRBs and lowest 
mean velocities for SNe Ic is present at almost all epochs. 

• The biggest difference in mean velocities for the dif¬ 
ferent sub-types of SNe Ic is at the earliest phases (at 
tv max = —10 days) and the difference decreases with 
time. Indeed, after tvmax=^0 days the Doppler absorp¬ 
tion velocities of all three subgroups become similar and 
are consistent within the standard deviation of the data. 
In general, SNe Ic have the smallest dispersion (exclud¬ 
ing PTF 12gzk), especially at later phases, while SN- 
GRBs have the largest dispersion - however since that 
dispersion is mostly driven by SN 2010bh/GRB 100316D 
, when discounting that SN-GRB the group with the 
largest dispersion (driven by the behavior of many SNe) 
is the group of SNe Ic-bl without observed GRBs. This 

We note that our conclusions hold even if our new method of 
estimating absorption velocities for SN Ic-bl via MCMC template 
fitting (Appendix [A|| is not correct, since the velocities measured 
via the classical identification methods are even higher than those 
from template fitting. 
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Phase since V-band maximum [days] 


Phase since V-band maximum [days] 


Fig. 5.— Fe II A5169 absorption velocities of different kinds of SNe, color-coded by SN type: SNe Ic-bl with observed GRBs, including 
the radio-relativistic SN2009bb, are denoted by filled red circles, SNe Ic-bl (without observed GRBs) by open blue circles, SNe Ic by filled 
grey circles. In addition we also show the velocity evolution of PTF 12gzk, which does not exhibit broad-lined spectra, but nevertheless 
shows very high absorption velocities, comparable to those of SNe Ic-bl. Left: measured values for individual SNe. Right: the rolling 
weighted mean of the individual SN measurements for their SN subtype, not including SN 2010bh/GRB 100316D nor PTF i2gzk, in rolling 
window sizes of 5 days (for phases with tymax < 30 days) or of 10 days (for phases with tyj^ax > 30 days), in steps of 1 day. The gaps 
are due to the fact that no weighted average was computed for fewer than three data points in the specific windows. The error bars for 
the weighted mean values are the weighted standard deviations of the data. The x-axis error bars on the phases indicate the average of 
the phases for the velocity measurement averaged in that windows. For reference, the phases H-IO, 0 and —10 days with respect to V-band 
maximum are marked. SN Ic-bl with GRBs have systematically higher velocities, followed by SN Ic-bl without observed GRBs, which have 
lower velocities, while SN Ic have systematically the lowest velocities - this conclusion holds even when excluding SN2010bh/GRB100316D 
(see right plot), which appears to exhibit the highest velocities by more than 150% than any other SN Ic-bl with a GRB. SN Ic-bl 2010ay 
is discussed in the text. 


TABLE 4 

Weighted mean absorption velocities 


SN Type 

II A5169)^ [km s 
at tvmaa;—0 days 

II A5169) [km s 
over all epochs 

SN Ic*^ 

-9,700 ± 1,500 

-8,900 ± 3,300 

SN Ic-bl (all) 

-21,100 ± 8,200 

-18,000 ± 8,000 

SN Ic-bl (no GRBs'=) 

-16,700 ± 2,700 

-13,300 ± 4,600 

SN Ic-bl -1- GRBs 

-23,800 ± 9,500 

-21,000 ± 8,000 


When excluding SN 2010bh/GRB 100316D 
SN Ic-bl (all) -18,500 ± 3,300 -16,200 ± 5,200 

SN Ic-bl -h GRBs -19,200 ± 2,900 -17,000 ± 4,000 


^ The weighted means were taking over a 5 day rolling windows centered at tv— 0, 
i.e. V-band maximum. The error bars refer to the weighted standard deviation of the 
data. For multiple velocity measurements for the same SN across the 5-day window, the 
weighted average for that SN was taken, as to not bias the mean SN type values by a 
few SNe that have many data points. Note that PTF 12gzk is not included in the mean 
values for SNe Ic. 

Excluding PTF12gzk. 

^ While no GRB emission was observed for these SNe Ic-bl, the necessary data to rule 
out an off-axis and/or low-luminosity GRBs are present for some, but not for all SNe 
Ic-bl in this sample. 


is the case at many phases, especially at 0< tvmax < 30 
days. 

• SN Ic-bl 2010ay shows high Fe II A5I69 absorption 
velocities, much higher th an those for SNe Ic-bl without 
observed GRBs (see also Sanders et al. 2012b). While 
no gamma-rays were detected trom SN 20l0ay, the 
gamma-ray and radio upper limits cannot rule out the 
possibility that it hoste d a common and low- luminosity 
GRB like GRB0602I8 (jSanders et al.||20I2b|. In fact, 
many of the properties of SIN 2010ay are similar to those 


of SN-GRBs: its high absorption velocity (sho wn here), 
its high peak luminosity, and its low metallicity (Sanders 


iiS _ 

et al.|2012b[ ) . Thus, we speculate that SN Ic-bl 2010ay is 


a very good candidate for a low-luminosity GRB, simi¬ 
lar to SN 2006aj/GRB 060218, seen either on- or off-axis. 


4.5. Line Widths in SNe Ic-bl with and without GRBs 

We present a novel approach to characterize and quan¬ 
tify the “broadness” of lines in SN Ic-bl spectra. Besides 
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the Doppler velocity, v, our template fitting procedure 
for SNe Ic-bl (Appendix also outputs the width of 
the Gaussian,cr, with whicn the appropriate SN Ic tem¬ 
plate has to be convolved to match the observed SN Ic-bl 
spectrum. We converted a to FHWM (which is straight¬ 
forward for a Gaussian) and call it the “FWHM convolu¬ 
tion” velocity. It thus indicates the width of the SN Ic-bl 
spectral line with respect to the SN Ic template at the 
^same phase. While the absolute width decreases as a 
function of time for the same SN (similar to the absorp¬ 
tion velocities decreasing in time), we find that there is no 
decrease in the convolution FWHM as a function of time, 
since these values are with respect to a SN Ic template, 
whose width decreases with time. Thus, since there is 
no time evolution, we plot all convolution FWHM values 
of SNe Ic-bl with and without GRBs in Fig. in terms 
of cumulative fractions, and also plot the (asymmetric) 
error bars as horizontal lines. We report median values 
in Table [H 



Convolution FWHM velocity [10^ km s 


Fig. 6.— Cumulative fraction (solid lines) of Fe II A5169 convo¬ 
lution FWHMs (based on the sigma velocities, see Appendix [A|| , 
with respect to the appropriate SN Ic template, of SN Ic-bl wTch 
observed GRBs (in red) and without observed GRBs (in blue), and 
their error bars on the velocities. These convolution FWHM veloc¬ 
ities with which the template SN Ic had to be convolved to match 
the individual SNe Ic-bl spectra indicate the width of the spectral 
line with respect to the SN Ic template. SN Ic-bl with GRBs have 
higher cr velocities, i.e. larger Fe II A5169 line widths than SN 
Ic-bl without GRBs. A K-S test yields the small probability that 
0.8% that the FWHM convolution velocities of SN Ic-bl without 
GRBs are drawn from the same parent population as SN Ic-bl with 
GRBs. 


We find that SNe Ic-bl (both types) have Fe II A5I69 
widths that are around 8,800 km s~^ (median value) 
broader than those of SNe Ic (based on data from all 


TABLE 5 

Median convolution FWHM velocities 


SN Type 

'*^convolutionFWHM (Fe II A5169)^ 
[km s“^] 

Median values based on all data points*^ 

SN Ic-bl (no GRBsG 

8,200 (N=87) 

SN Ic-bl -I- GRBs 

9,300 (N=103) 

Median values based 

on the weighted mean for each SN 

SN Ic-bl (no GRBsG 

7,600 (N=10) 

SN Ic-bl -I- GRBs 

9,300 (N=9) 


^ These convolution FWHM velocities are with respect to the 
the width of the SN Ic template at the appropriate phase. 

^ May include multiple data points per SN. 

^ While no GRB emission was observed for these SNe Ic-bl, the 
necessary data to rule out an off-axis and/or low-luminosity 
GRBs are present for some, but not for all SNe Ic-bl in this 
sample. 


epochs). In addition, figure 1^ shows that SN-GRBs have 
consistently higher convolution Fe II A5169 FWHM ve¬ 
locities than SNe Ic-bl without GRBs by 1,000 — 2,000 
km s~^(see Table ^ - the same trend as for absorption 
velocities (Section T-4[ ). A Kolmogorov-Smirnov (KS) 
test indicates that there is a small chance, namely i.8% 
that the FWHM convolution velocities of SN Ic-bl with¬ 
out GRBs are drawn from the same parent population 
as SN Ic-bl with GRBs. As before there may be the con¬ 
cern that for this approach, namely including all avail¬ 
able data points for each for each SN, would bias the 
results towards the SNe with many data points. Thus, 
using the same procedure as in Section [44l we took the 
weighted average of all data points for each SN and then 
computed the median for all SNe of the same type. The 
results are unchanged, and the result of the KS test is 
0.3%. To get an estimate on how the velocity error bars 
(shown in Fig. ^ may influence the results, we compute 
the KS probability for 2 extreme scenarios using data at 
all phases: when the two velocity distribution are closest 
(i.e., adding the errors to values of the SN Ic-bl without 
GRBs and subtracting the errors from the SN-GRB val¬ 
ues), and when the two are furthest. For the scenario 
where they are closest, the KS test still obtains a low 
probability of 3%. Thus, even including error bars and 
the most extreme case, the line widths of SN-GRBs and 
those of SN Ic-bl without observed GRBs are unlikely to 
be drawn from the same parent population. Of course, 
for the scenario when they are furthest, the KS test ob¬ 
tains the 0% probability within our significance. 

In summary, we conclude that the Fe II A5I69 widths 
of SN-GRBs are higher than those of SNe Ic-bl without 
observed GRBs (by ~ 1,000—2,000 km s“^), in line with 
the visual trend for all features in the mea n sp ectra of 
SNe Ic-bl with and without GRBs (Sectionj3.1[ Fig. 
and similar to the tren d seen for Fe II A5169 absorption 
velocities (Section 4.41. We discuss the implications of 
this finding, together with our speculations, in the next 
section. 


4.6. Absorption Velocities vs. Line Width Velocities 

For all SNe Ic-bl in our sample, we plot the absorp¬ 
tion velocity against the line-width velocity based on the 
same ^ectrum, measured as close as possible to tvmax 
(Fig. [^. We note that for some SNe, the spectrum clos- 
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est to tvmax IS actually late, e.g., at tymax=15.7 days 
for SN 2010ay. While in many cases, high absorption 
velocities correlate with broad lines (the Spearman coef¬ 
ficient is p = —0.599 with a strong significance of 0.004, 
and the Pearson coefficient is r = —0.642), there are two 
SNe for which that is not the case (SNe PTFlOvgv and 
PTF12gzk). This is the first time such a correlation has 
been investigated explicitly and methodically for a large 
sample of SNe Ic-bl. 



Absorption velocity [10^ km s '] 


Fig. 7. — Comparison between Fe II A5169 convolution FWHMs 
and absorption velocities with respect to the appropriate SN Ic 
template, of SNe Ic-bl with observed GRBs (in red; in addition, 
filled circles: with LLGRBs and their names reported; open cir¬ 
cles: with HLGRBs) and without observed GRBs (in blue) for 
the spectrum that is closest to the respective SN date of V-band 
maximum. These convolution FWHM velocities with which the 
template SN Ic had to be convolved to match the individual SNe 
Ic-bl spectra indicate the width of the spectral line with respect 
to the SN Ic template. For reference we also mark the weighted 
average absorption velocity for SNe Ic (dashed line), as well as 
its standard deviation (dotted lines), and its convolution velocity 
of 0 di 1000 km s“^. While absorption velocity and width ve¬ 
locity are correlated for many SNe, there are also outliers, such as 
PTF 12gzk, which has narrow lines (it is a SN Ic), but high absorp¬ 
tion velocity, and PTF lOvgv, whose absorption velocity is lower 
than the average for SNe Ic, but still is a broad-lined SN Ic. We 
also plot in orange the values for the radio-relativistic SNe 2009bb 
and 2012ap. See text for discussion on any differences between 
LLGRBs and HLGRBs. 


As shown. PTF 12gzk possesses very high absorption 
velocitiesjlj well within those of SN Ic-bl with GRBs, 
but narrow lines (sim ilar to those of SNe Ic, see also 
Ben-Ami et al. |20I2 ). In terms of reasons for the mis- 


match between the two kinds of velocities for PTF 12gzk, 


Ben-Ami et al.| (2012) speculated that it could be due 
to aspherical explosion geometry, or high ejecta mass or 
very steep density gradient. PTFlOvgv, on the other 
hand, possesses broad lines, but has the lowest absorp¬ 
tion velocity amongst all SNe Ic-bl, lower by a factor 
of ^1.5 than the average value for SNe Ic-bl, and even 
lower th an the average for SNe Ic. This must be the rea¬ 
son why Corsi et al. (2011) classified this SN as a SN Ic, 


Radio observations suggest that the synchrotron radiation 
emitting layer, i.e. the fastest-moving ejecta, also possessed high 
veloc ities, with a mean velocity of 80,000 kms~^(Horesh et al. 


[Mi^. 


since their definition of SNe Ic-bl relies on the value of 
the “photospheric” velocity (in their case, inferred from 


spectral s ynthesis modeling based on the code from Maz- 


zali 2000), which is low for PTF lOvgv. However, the 


dehnition of a broad-lined SN Ic relies on the measured 
width of the line and should be adopted for classifying 
SNe. 

We highly encourage modeling the spectra of SNe Ic- 
bl with a detailed parameter study in terms of e.g., as- 
phericities, density gradients and ejecta masses, to re¬ 
produce the observed parameter space of all SNe Ic-bl 
and not to just fit a hand-crafted model spectrum to an 
individual SN. 

In addition, we also search for differences in SN Ic- 
bl spectra for those connected with Low-Luminosity 
GRBs vs those connected with High-Luminosity GRBs, 
as some theories suggest different physical mechanisms 
for LLGRBs and HLGRBs. Figures and sug¬ 
gest no obvious trend between the SN spectralprop- 
erty (specihcally, Vabs) and GRB luminosity. While 
SNe Ic-bl connected with LLGRBs appear to show a 
larger range in velocities (—16,000 < Vabs < —40,0000 
km s~^ at V-band maximum) than SN Ic-bl connected 
with HLGRBs (—19,000 < Vabs < —27,0000 km s“^ 
at P-band maximum), this is only due to one object, 
SN 2010bh/GRB I00316D, which has velocities that are 
much larger, up to 150%, compared to the rest of SN- 
GRB. In other words, for velocities at P-band max¬ 
imum, the standard deviation of the HLGRB sample 
is ~ 3300 km s~^, while that of the LLGRB with SN 
20I0bh/GRB I003I6D is ~10,000 km s~^, however with¬ 
out it, the standard deviation drops to ^2500 km s“^, so 
becomes similar to that of the HLGRB sample. Given 
the small data set of a total of four objects, it is hard 
to determine whether SN 2010bh/GRB I003I6D is an 
outlier, and thus more data are needed. 

We further discuss the i mpl ications of this work in Sec¬ 
tion]^ and specifically in|5.3[ 


5. DISCUSSION: CONSTRAINTS ON PROGENITORS AND 
EXPLOSION MECHANISMS 

Although only SNe Ic-bl have been connected with 
GRBs, the majority of SNe Ic-bl show no detectable 
gamma-ray emission, and their relationship to normal 
SN Ic is unknown. Thus, one of the key outstanding 
questions is what conditions lead to each kind of explo¬ 
sion in which kinds of massive stripped stars. We have 
explored this question by investigating what sets apart 
the spectra of SNe Ic-bl with GRBs from those without 
GRBs, and how SNe Ic-bl spectra are related to those 
of SNe Ic. We have uncovered a number of trends while 
analyzing our datasets. Here, we discuss those that are 
specifically related to constraining progenitors and ex¬ 
plosion mechanisms. 


5.1. Helium in Proqenitors of SN Ic-bl and of 
SN-GRBs? 

In Section |3.2| we addressed the question whether He 
lines could be smeared out in the spectra of SNe Ic-bl 
such that they would not be detectable, despite He be¬ 
ing present in the ejecta. The lack of He-detection in 
optical spectra of SNe Ic-bl and GRB-SNe is puzzling, 
since many of the progenitor models, especially those of 
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SN-GRBs, are supposed to be He-stars. Here, we sim¬ 
ulated a “broad-lined SN Ib”, i.e., a broadened SN Ib 
with He lines. While convolving a SN Ic mean spectrum 
with a Gaussian profile yields a spectrum that is very 
similar to a SN Ic-bl mean spectrum, a similarly con¬ 
volved SN Ib mean spectrum (effectively a “SN Ib-bl”) 
is different from the SN Ic-bl spectrum by more than 2a 
over the wavelength ranges of interest. We tested this 
for the most extreme cases of smearing (i.e., for large 
convolution velocities). Thus, we conclude that the ab¬ 
sence of clear He lines in optical spectra of SNe Ic-bl and 
of SN-GRBs is not due to the smearing out by the high 
velocities present in SNe Ic-bl. This implies that the pro¬ 
genitor stars of SN-GRBs may be truly free of He, in ad¬ 
dition to being H-free, putting strong constraints on the 
stellar evolutionary paths needed to produce such GRB 


licities (e.g, Modiaz et al. 

2008 

Sanders et al. 

2012a 

Graham & I'ruchter 2013), including the potential need 


and therefore burn it (Frey et al. 2013). 


Since standard mechanisms in stellar evolutionary 
models for single stars find it very difficult to remove 
so much He such that no or very little He (less than 
0.2 Mq ) remains (e.g., Yoon et al. 2010 Yoon 2015), 
especially at the lower nietallicities, at which SiN-GRBs 
are generally found, one of the only remaining paths 
may be rapid mass loss after core helium exhaustion 
due to interaction with a compact companion in a 
very short-peri od binary system (S-G. Yoon, private 


communication, 


Pols & Dewi 2002 Ivanova et^jdj 2003| 


in potentially young dense star cluster (van den Heuvel 


fc Portegies Zwart 2013) or very miassive He-stars 


(| Yoon et al. 2010 see their Fig. 10, but only for solar 


metallicityj. Xnother viable alternative is explosive 
common-envelope ejection in a slow merger of a massive 
primary with a low mass companion (?). Thus, our 
He-constraint supports short-period binary models for 
SN-GRBs. This suggestion is in line with recent obser¬ 
vational work on host gala x ies of SNe Ic-bl, SN-GRBs 


and GRBs by Kelly et al. (2014), who interpret their 
preference for over-dense galaxies as indicating that 
they preferentially come from tight binary systems, 
which may be produced more readily in such over-dense 
conditions. There is yet another mechanism to “hide” 
He in SN Ic spectra: if the amount of ^®Ni mixing is 


the optical spectra 

(Lucy 

1999 

Hachinger et al. 

2012 

Dessart et al. 

2012) 

. However, 

based on our thorough 


an alysis ot the larg est data set of SNe Ib and SNe Ic 


Liu et al. (2015), we think that there is no hidden 


He in normal sNe Ic , since the SN observables are 


inconsistent with the predictions of Dessart et al. (2012), 
who had modeled the effect of insutticient mixing in SNe 
Ic spectra that may have lead to hiding He. Similar 
conclusions were reached by [Cano et al.| 2014| based on 
modeling of SN Ib 1999dn arid SN Ic 2007gr. 

Nevertheless, in order to fully understand the progeni¬ 
tors of SNe Ic-bl, we highly encourage modeling the spec¬ 
tra of SNe Ic-bl in detail to understand the physical con¬ 
ditions and to do a detailed parameter study, which has 
not been done before, using sophi sticated spectral syn - 
thesis models such as GFMGEN (Dessart et al. 2011), 


PHOENIX (Baron et al. 1995), and SEDONA (Kasen 
et al. 2006).^ In this context, we note that there are 
discrepancies between the simulated SN Ic-bl spectrum 
(broadened SN Ic mean spectrum) and the real SN Ic-bl 
spectrum for two different regions: a) in the blue re¬ 
gion (3500-4400 A, Fig. |^, where the forest of Fe lines 
lie, which may have implications for understanding line 
formation in highly-expanding atmospheres; and b) the 
region around 6300A (with absorption commonly associ¬ 
ated with Si II A6355). In addition, we find that high- 
absorption-velocity SNe may not have broad lines (e.g., 
PTF12gzk), and conversely, that broad-lined SNe Ic can 
have low absorption velocities (e.g., PTFlOvgv). Given 
that these two SN mavericks were found by PTF, we are 
looking forward to more interesting SNe discovered via 
untargetted and innovative surveys, such as LSST. In the 
meantime, until a unifying picture of spectral line for¬ 
mation in SNe Ic and SNe Ic-bl is developed and tested, 
we suggest that authors, when reporting velocities, are 
clear about how the number was derived (whether from 
modeling or observations, and which line was used) and 
whether it is the absorption velocity or the width of the 
line they are reporting, and to use the same method when 
comparing different SNe. 

5.2. Progenitors and Explosion Conditions of SNe Ic-bl 
with GRBs vs. SNe Ic-bl without observed GRBs 

We have found that SNe Ic-bl with obse rved jets have 
quantifiably broader lines (Sections |3.3| a nd 4.6) and 
higher absorption velocities (Section |4.4), and there- 
fore higher ejecta velocities than those without observed 
GRBs. One explanation is that the latter either had 
lower-energy or shorter-duration jets that were stifled 
and could not break out of the star, but neverthe¬ 
less imparted some of their energ y to the stellar enve - 
lope and therefore accelerated it (Lazzati et al. |2012 ), 
though to lower velocities than for SiN-GRBs. In this 
scenario, radio-relativistic SNe such as SNe 2009bb an d 
2012ap (jSoderberg et al.||2010| jMargutti et al.||2014|), 
which are claimed to be engine-driven, but with no 
GRB detection, are explosions where the jet may have 
lasted for a time that is shorter than for the SN-GRB 


(Margutti et al. 2014), but longer than for a SN Ic-bl 
without a GRB. Another reason why the jet could have 
been choked is because the stellar envelope, which be¬ 
came SN ejecta when the star exploded, had a larger 
mass in SNe Ic-bl without GRBs, and thus gave rise to 
lower velocities, as observed, for the same kinetic en¬ 
ergy as in SNe Ic-bl with GRBs. If both kinds of SNe 
Ic-bl had the same kinetic energy, the ejecta mass of 
SNe Ic-bl without observed GRBs would be on aver¬ 
age < VabsiSNGRB) >^ / < VabsiSNIc — bl) >^= 2.6 
times higher compared to the ejecta mass of SN-GRBs 
(or twice when excluding SN 2010bh/GRB 100316D). 
This extra mass may lead to the stifling of the jet. Con¬ 
versely, if both SN-GRB and SNe Ic-bl had the same 
ejecta masses, but released different amounts of kinetic 
energy, then SN-GRBs would be releasing 2.6 times more 
energy than SNe Ic-bl without GRBs (or twice if ex¬ 
cluding SN 2010bh/GRB 10 0316D). At face value, spec- 
tral synthesis modeling by Mazzali et al.| (2013) and 
collaborators suggests that SiN 19y8bw/GRB 980425, 
SN 2003dh/GRB 030329, SN 20031w/GRB03I203 have 
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higher ejecta masses and higher explosion energies than 
SNe Ic-bl without GRBs, based on their small sample 
(though keeping in mind all the caveats about geometri¬ 
cal assumptions and sim plifications). This result is in 
line with calculations by Cano (2013) on a larger set 
of SNe Ic-bl and SN-GRBs, though th ere are a nu m- 
ber of assumptions and simplifications in |Cano||2013[ in¬ 
cluding using Fe- and Si-based velocities interchangeably, 
even though they give systematically different values as 
shown above. Thus, we conclude that it is reasonable 
that SNe Ic-bl without observed GRBs, which appear to 
have intrinsically lower kinetic-energies than SN-GRBs, 
may have had also lower energy jets that were stifled in 
the passage through the star. In addition, it is possible 
that both phenomena (successful jet and higher veloci¬ 
ties in the SN ejecta) are caused by th e same ingredient 
such as high angu lar momentum (e.g. Fuller et ar]|2015 


Mosta et al.||201^. There is an assumption in sortie pa- 


pers that the GRB jet inside the star has on average a 
small energy of ^ 10®^ ergs (e.g., Mazzali et al. 2014). 
This is a misconception, since lO'^^ ergs is the average 
gamma-ray radiation energy of the jet after it breaks out 
of the star, which can be lower due to losses and ineffi¬ 
ciencies in producing radiation. Indeed, most simulations 
include much stronger jets that travel throug h the star. 
almost always with energies of 10®^ ergs (e.g., Lazzati et 
‘M^lDuffell fc MacFadyen||2015| ). 


al. 


'hus, it is reasonable to conclude that SNe Ic-bl with¬ 
out observed GRBs, which appear to have intrinsically 
lower kinetic energies than SN-GRBs, may also have had 
lower-energy jets that were stifled in the passage through 
the star, potentially more affected by the kink instability 
( Bromberg fc Tchekhovskoy||2016 ). 


SN-GRB central engines (e 

g., Usov 

1992 Thom 

Dson et 

al. 

2004 

Bucciantini et al. 

12007" 

M 

etzger et al. 

20111), 

whose rotational energy reservoir of 10°^ ergs (Usov|I992 


but see more recently 

Metzger et al. 

2015) appears 

be consistent with the 

tliffh average 

UKB-yiN energies 

^ 2 X 10^^ ergs ( 

Mazzali et al.|20I4 

Cano et al.|20I5 ) 


to be maximally efficient in extracting the rotational en¬ 
ergy of the rapidly rotating pulsar, which is somewhat 
hard to imagine given the usual inefficiencies in astro- 
physical systems. 

Another potential explanation is that the lower ve¬ 
locities of SNe Ic-bl without observed GRBs is due to 
viewing angle effects, where SNe Ic-bl have no observed 
GRB because their beamed GRB jets were not point¬ 
ing along our line of sight; such SNe harboring off-axis 
GRBs would presumably also have lines that are less 
broadened, since we would not be looking at the SN 
ejecta that got most accelerated by the passage of the 
jet along the line-of-sight. We think that it is unlikely 
that all nearby SNe Ic-bl without observed GRBs in our 
sample are due to off-axis GRBs, for the following rea¬ 
sons: more than half of the nearby (z < 0.3) SN-GRBs 
have very large opening angles of u p to 180 deg (see 
e.g., Table I in Mazzali et al. 2014), such that highly 


collimated GRBs are probably rare in the nearby uni- 

Ho wever, note that the KE was estimated in various dif ferent 
wavs i n llMa zzali et al.|2014[l and that the SN KE values in|Cano| 
|et al.| | | 20 i&| are not corrected for asymmetry. 


verse (e.g., Guetta & Della Valle 2007 Baerwald et al. 
20I5|). In addition, other lines of evideiice suggest that 


the majority of SNe Ic-bl are not off-axis GRBs: deep 
radio observations have not dete cted off-axis GRBs in 
any SN Ic-bl or any Stripped SN ( Soderberg et al.|[2006 


|20I0[ jCorsi et ^|20II[ |20I2[ R. Margutti, private com- 
munication 2015) — and in the cases where radio emis¬ 
sion was detected, it was consistent with a model of the 
pure hydrodynamic explosion, except for the t wo cases of 
engine-driven SNe Ic-bl 2009bb and 20I2ap (Soderberg 

metallicities 


et al. 


bl wf 


20101 IMargutti et al.||20I4|). Furthermore, SNe I 
fh gRBs are at systematically lower metalliciti 


than SN Ic-bl without observed GRBs (Modjaz et al. 
2008 Modjaz 2012 Sanders et al.||20I2a ), which cannot 


be explained by viewing angle effects. 'I’hus, it is unlikely 
that all nearby SNe Ic-bl without observed GRBs in our 
sample are due to off-axis GRBs, and thus we favor the 
first hypothesis, namely that SNe Ic-bl without observed 
GRBs had lower-energy or shorter-duration jets, as also 
supported by the computed kinetic energies of a few SNe 
Ic-bl, compared to those of a few SN-GRBs. 

In either case, both scenarios support the hypothesis 
that the engines of SNe Ic-bl (with and without GRBs) 
are fundamentally different from those of SNe Ic, and 
may come from different physical progenitors than SNe 
Ic. Indeed, other lines of evidence point in the same 
direction: SNe Ic-bl and GRBs prefer over-dense host 
galaxy environme nts , compared to other Stripped SNe, 
including SNe Ic ( Kelly et al.|20I4 ), and SNe Ic-bl (with 


ities than those of SNe Ic 

Modiaz et al. 

2008 

ModiazI 

2012 Sanders et al.||20I2a 

Kelly & Kirshner|201 



5.3. Are SNe le-bl connected with LLGRBs different 
from those conneeted with HLGRB? 

Starting with the first case of the SN-GRB connec¬ 
tion, namely SN I998bw/GRB 980425, there have been 
so far four SN-LLGRBs observed, which as a class. 


(e-So 

Cobb et 

al. 

2006 Soderberg et al. 

2006 

Guetta 

& Delia Valle 

20 

J7). In addition to their low lumi- 


more isotropic and have relativel y soft high-energy spec¬ 


tra compared to HLGRB s (e.g., Campana et al. 


2006 


Hjorth| 2013 Nakar|[2M5 and references therein), 'there 


have been suggestions that the high-energy emission of 
the GRB in these two groupmay come from different 


mechanisms 


2011 




Waxman et al. 


Nakar & Sari 


: breakout of 


i 


2007 Bromberg et al. 


T5), namely that LL- 


GRB are from shock breakout of the collapsar jet from a 
low-density envelope or wind, and thus are more spheri¬ 
cal because of that. 

In figure we plot the measured SN absorption veloc¬ 
ity as a function of isotropic gamma-ray luminosity (left) 
and of gamma-ray energy (right ), showing both isotropic 
and bea ming-corrected energ ies (taken from the compi¬ 
lation in MazzaIj_gt^ay|20I4l see references therein, and 
also see Ryan et al. 2015). The velocities are based on 


the spectrum within tvmax = i2 days. As can be seen, 
when corrected for beaming angle, the intrinsic energies 

We note that there are also some GRBs whose 7 -ray lumi¬ 
nosity appear to be intermediate (e.g., SN 2012bz/GRB 120422A, 
Schulze et al. 2014). 





















































































































































19 


_ -45 
c 

2 -35 
-30 

U 

o 

^ -25 
a 

■a -20 

& 

o 


-15 

-10 



• SN Ic-bl+LLGRB 
O SN Ic-bl + non-LLGRB 










48 49 50 51 52 

log Lyiso [erg] 


CJ 

"oj 

> 

C 

o 


< 


-45 

-40 

-35 

-30 

-25 

-20 

-15 

-10 


O 

□ 


SN Ic-bl+LLGRB 
SN Ic-bl + non-LLGRB: 

p ; 

yconecfed 



48 49 50 51 52 53 54 

logEy.,,, [erg] 


Fig. 8. — Comparison between SN spectral properties and the high-energy properties of their accompanying GRBs: absorption velocity 
(of Fe II A5169 at V-band maximum) as a function of isotropic gamma-ray luminosity (left) and of gamma-ray energy {rights showing both 
isotropic and beaming-corrected energies, connected by a dashed line for beamed GRBs). The velocities are based on the SN spectrum 
within tvmax = days. No strong correlation between SN spectral properties and GRB properties are seen (see text for more discussion). 


are similar between LLGRBs and HLGRBs and the only 
difference appears to be collimation, with LLGRBs being 
almost spherical. There is no clear correlation between 
the SN absorption velocities and the GRB isotropic en¬ 
ergies and luminosities (for both cases, the Spearman 
coefficient is p= 0.2 with a low significance of 0.6), nor 
beaming-corrected (the Spearman coefficient is p= 
0.4 with a low signihcance of 0.3). A similar analysis 
correlating the SN peak optical lumino sity with its G RB 
luminosity did not find a trend either ( Hjorth|[2013 ). 

6 . CONCLUSIONS & FUTURE DIRECTIONS 

We have presented the first systematic investiga¬ 
tion of spectral properties of explosions from massive 
stripped stars, which are sometimes accompanied by 
long-duration GRBs, in order to probe their explosion 
conditions and their progenitors. To that aim, we ana¬ 
lyzed a total of 407 spectra of 11 SNe Ic-bl with GRBs, 
10 SNe Ic-bl without observed GRBs, and 17 SNe Ic. 
The data consist of published spectra of individual SNe 
from the literature as well as the densely time-sampled 
and homogeneous sp ectral data set of stripped SNe from 
Modjaz et al.|(|2014). 

In order to quantify the diversity of the spectra of a 
specific sub-type, we constructed average spectra of nor¬ 
mal SNe Ic, SNe Ic-bl with GRBs and SNe Ic-bl with¬ 
out GRBs, along with standard deviation and maximum 
deviation contours. We furthermore explore the ejecta 
velocities present in different kinds of SNe Ic using two 
different methods: 1) measuring Fe II A5169 line veloci¬ 
ties from their Doppler-shifts and 2) exploring the widths 
of the lines. We use a lot of care to derive robust error 
bars for all our measurements by constructing spectral 
uncertainty arrays for all our spectra from the SN spec¬ 
tra themselves and by propagating them into our final 
measurements, either via Monte Carlo resampling or by 
formal fitting using MCMC. In addition, we have devel¬ 
oped a novel method for analyzing the blended features 
of SNe Ic-bl, where in particular the three lines of Fe II 
(Fe II A4924, Fe II A5018, and Fe II A5169) are blended 
into one broad feature in SNe Ic-bl. This novel method 


uses our prior constructed mean SN Ic spectra as tem¬ 
plates for the appropriate phase and hnds a best fit to 
the SN Ic-bl spectra, using MCMC to explore the full 
parameter space and to derive error bars. 

The main results that we have obtained can be sum¬ 
marized as follows: 

• SN Ic 19941 is not a typical SN Ic, in contrast to com¬ 
mon usage, while the spectra of SN 1998bw/GRB 980425 
are representative of mean spectra of SNe Ic-bl. While 
the spectra of SN 1998bw/GRB 980425 lie close to the 
mean spectra of SNe Ic-bl with GRBs and well within 
one standard deviation, those of SN 19941 are close to 
the maximum deviation of the mean spectra of SNe Ic 
(especially at tvmax=^0 days) at many wavelengths. 

• For SN Ic 19941 in particular, and for all SNe in 
general, the velocity obtained depends on which method 
was used to measure photospheric velocity, as differ¬ 
ent lines (Fe II A5169 vs. a line commonly identified 
as Si II A6355) and different spectral synthesis codes 
(SYNOW vs. Mazzali’s NLTE code) yield systematically 
different values. Furthermore, it is crucial to distinguish 
between two kinds of velocities (absorption velocity vs. 
line-width velocity) since they are not always correlated 
for the same SNe (e.g., PTF 12gzk and PTF lOvgv), 
and we encourage authors to be clear about which kind 
of velocity they are reporting, as well as to use the same 
method when comparing different SNe. We encourage fu¬ 
ture modeling of the spectra of SNe Ic-bl with a detailed 
parameter study to reproduce the observed parameter 
space of SNe Ic and Ic-bl, including asphericity in the 
ejecta and viewing angle effects. 

• SNe Ic-bl (with and without GRBs) have Fe II A5169 
widths that are ^9,000 km s“^ (median value) broader 
than those of SNe Ic. They also have average 
Fe II A5169 absorption velocities that are larger than 
those of SNe Ic by ^9,000 km s~^ when excluding SN 
2010bh/GRB 100316D and by ^13,000 km s“^ when in¬ 
cluding SN 2010bh/GRB 100316D. The latter values are 
for the epoch of R-band maximum. 

• Somewhat surprisingly, we find that SNe Ic-bl with 
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GRBs, on average, have higher ejecta velocities (by 
-6,000 km if excluding SN 2010bh/GRB 100316D), 
as well as wider line widths (by — 2,000 km s~^) 
than SNe Ic-bl without observed GRBs, as traced by 
Fe II A5169. That is, the SN Ic-bl “knows” whether 
it harbored a GRB beamed towards us. We can inter¬ 
pret these observations in two ways: the first is that 
those SNe Ic-bl without observed GRBs had lower-energy 
or shorter-duration jets that were stifled and could not 
break out of the star, but nevertheless imparted some 
of their energy to the stellar envelope and therefore ac¬ 
celerated it less than SNe with high-energy, visible jets. 
The other possibility is that SNe Ic-bl without observed 
GRBs did harbor GRBs that are highly beamed, but 
they were all beamed away from us, i.e., they were all 
off-axis GRBs. We disfavor the latter possibility since 
there is strong evidence that a number of the SNe Ic-bl 
in our sample did not harbor off-axis GRBs, and since 
many SN-GRBs are not highly beamed in general. In 
order to verify our hypothesis we suggest that the geom¬ 
etry of the ejecta of SNe Ic-bl without GRBs should be 
probed (e.g., via polarization studies or late-time spec¬ 
tra) to detect any other signs of a stifled jet inside SNe 
Ic-bl without observed GRBs. 

• We search for differences in SN Ic-bl spectra for those 
connected with Low-Luminosity GRBs vs. those with 
High-Luminosity GRBs, as some theories suggest differ¬ 
ent physical mechanisms. While the sample sizes are still 
very small (with each group having —4-5 objects), we find 
no clear correlations between SN spectral properties and 
GRB energies, nor GRB luminosities. 

• The absence of clear He lines in optical spectra of 
SNe Ic-bl cannot be due to the smearing out of the He¬ 
lium lines by the high velocities present. While the lack 
of ®®Ni mixing may in principle be responsible for ’’hid¬ 
ing' 


He l ines, we show in our companion paper |Liu et| 
al. (2015) that the lack of ®®Ni mixing is most likely not 


Hachinger et al.||2012 

Cano et al. 

2014). These findings 

imply that the progenitor stars o 

r SJNe Ic-bl (with and 


without GRBs) are probably He-free, in addition to be¬ 
ing H-free, which puts strong constraints on the stellar 
evolutionary paths needed to produce H-free and He-free 
progenitor stars at low metallicity. In particular, our He 
constraint supports short-period binary mo dels for SNe 


Ic-bl. in line with other recent works (e.g., Kelly et al. 


2014). In light of this, we encourage more studies of star 


formation conditions and of the shape of the initial mass 
function in galaxies similar to those host galaxies of SNe 
Ic-bl and GRBs, namely over-dense, low-mass, and low- 
metallicity galaxies, as well as theoretical investigations 
of removing both H and He in the potential binary sys¬ 
tems of GRB progenitors. 


• Overall, we encourage more observations of SNe Ic-bl 
without observed GRBs - while we have similar numbers 
of SNe Ic-bl with and without GRBs in this work, the 
number of spectra for SN-GRBs is a factor of —3 larger 
than those of SNe Ic-bl without observed GRBs. This 
is also in contrast to the fact that many more SN Ic- 
bl are known to be without observed GRBs than with 
observed GRBs. In addition, deeper radio and X-ray 
limits are needed for SNe Ic-bl to test for the presence of 
low-luminosity GRBs and low-energy outflows. Recently, 
there have been some suggestions that SN-GRBs may 
have standardiz able light curves (Gano et al. 2014 O 


& Hjorth 2014), which may have implications on any 
relationships between ®®Ni production (which sets the 
luminosity) in the different GRB engine scenarios and 
the physics of the SN ejecta, including its geometry and 
the observer’s viewing angle. 

All data products and codes are open access and 
open source and are available on our SNYU github 
website^^ including the code and MGMG wrapper for 
measuring velocities in SNe Ic-bl, and the additional 


as in|Liu & Modjaz 

(2014 

2(J14f are available 

on tr 
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APPENDIX 

A. MEASURING LINE VELOCITIES VIA THE TEMPLATE-FITTING METHOD IN SNE IC-BL 


The blueshift and width of specific spectral absorption features, such as Fe II A 5169, are used to track the photo- 
spheric velocity of SNe. For SN Ic spectra, it is straight-forward to identify absorption due to Fe II A5169, while Fe II 
AA 4924, 5018 lines are usually blended into a single feature (see Fig. |^. For SNe Ic-bl, all three Fe II lines are heavily 
blended due to the high expansion velocities in the ejecta (see Fig.[^and[l0|. Thus, it is difficult to accurately identify 
Fe II A 5 169 and measure its absorption velocity. There have been several attempts in the literature to solve this 
problem. Lyman et al. (2014), for example, used Si II A6355 instead, whenever they could not identify Fe II A5169. 

I I A5 169 absorption velocity is systematically higher than the Si II A6355 measured using 


Lyman et al. (2014) 
However, we hnd that the he 


the same method (see Section 4.4 and Bianco et al, in prep). Thus, the two velocities cannot be used interchangeably 
to track photospheric velocity m order to calculate physical parameters for SNe. What makes it worse is that Si II 
A6355 is sometimes blended with the Na I D lines in t o one absorption feature (e.g., PTF lOqts). Thus, it is at times 
difficult to identify Si II A6355 as well. Schulze et al. (|2014) analyzed spectra of SN -GRBs (many of the same spectra 
as in our SN-GRBs sample) and attributed the blended and broad feature du e to the three Fe II l ines to only Fe II 

up to —50,000 km s“^ 


Pignata et al. 


(2011) used SYNOW, a 


A5169. Thus, they derive systematically high velocities, 
spectral synthesis code, to obtain velocities of Fe II AAA4924, 5018, 5169 for Sf^ 2U09bb. Uitf'erent sp ectra l synthesis 
codes, however, yield different velocities for the same SN using different methods, as shown in Section [4)31 
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Fig. 9.— SN spectra of a SN Ic and a SN Ic-bl showing the wavelength region around the Fe II lines. 


We developed a data-driven and consistent method to identify Fe II A5169 and measure its velocity in the blended 
spectra of SNe Ic-bl: we used a template fitting approach in which we found the best fit between a convolved and 
blueshifted SN Ic template and the SN Ic-bl spectrum under consideration, at simil ar phases, for the wavelength region 
of the feature under consideration. This is similar to what we did in Section |3.2| for the whole spectrum, but now it 
is a formal fit, which we implement via Monte Carlo Markov Chain (MCMC) methods. The fits that this procedure 
yields are generally good, as indicated by the ~ unity values for the reduced values. 

Our SN Ic template spectra are the mean SN Ic spectra of Section but with a different time spacing: no w w e 
construct them for every two-day-phase interval (as opposed to five-day-intervals for the mean spectra of Section 3.1), 
in the range between t = —10 and t = -1-72 days with respect to maximum light (i.e, t = —10, —8, —6, —4, etc. days). 
Each mean spectrum includes SN Ic spectra in 5-day bins centered around the middle, i.e. ’’target” phase. For example: 
in order to measure the photospheric velocity for a SN Ic-bl spectrum at phase t =0 days, we use the mean SN Ic 
spectrum constructed using SN Ic spectra between phases —2.5 and 2.5 days and centered at phase t — 0 days. Thus, 
we note that template spectra at adjacent phases may not be independent from each other, but this aspect should not 
influence our results. For measuring the velocity of a SN Ic-bl spectrum at phase t = —2 days, we use the mean SN 
Ic spectrum centered at t = —2 days (constructed from SN Ic spectra between —4 and 0 days) for template fitting. 
Each mean spectrum includes no more than one spectrum per SN Ic, in order to avoid biasing the mean spectrum 
towards better-observed SNe. If two or more spectra of the same SN Ic are available within the phase range (within 
±2.5 days of the target phase), we choose the one nearest to the target phase. If there is still more than one spectrum 
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Rest wavelength (A) Rest wavelength (A) 

Fig. 10.— Example for convolution fitting of the Fe II feature for the spectrum of a SN Ic-bl (here: PTFlOqts) at t=0 days. In both 
plots, the SN Ic-bl spectrum is in black with its error spectrum overlaid in grey. Left: The SN Ic template spectrum (before the fit) is 
overplotted in red, and shows a “W” feature between 4600 A and 5300 A due to the Fell triplet: the red absorption trough is due to Fe II 
A5169, while the blue one is a blend of the other two Fe II lines. In the spectrum of the SN Ic-bl PTFlOqts, the three Fe II lines are blended 
into one feature. Right: The fitted SN Ic template spectrum is overplotted in red, with the 50th percentile confidence values and their error 
bars (corresponding to 50*^ — 16*^ and 84^^ — 50^^ percentiles respectively) determined by Markov-Chain Monte Carlo (MCMC) methods. 
The SN Ic template has been dopplershifted by the 50th percentile value of velocity parameter v (= 11,260(|(280 broadened by 

the Gaussian kernel with the 50th percentile value of sigma parameter cr (= 4460^32g km s“^), and stretched along the y-axis with the 
best value of amplitude parameter a (= 1.4 ^q'3 ). The thick red line denotes the spectrum that is used for the fit with the best value of the 
wavelength-range parameter Am (= 2.5^2 q t>y which the initial wavelength range (between the two maxima) gets subtracted. Note 
that the final reported value (and its uncertainty) of the Doppler velocity parameter v includes the contribution of the SN Ic template itself 
for each phase at which the fit is performed. Also shown on the plot is the reduced value for this fit. 


that satisfies the above condition, we use the one that covers the largest optical wavelength range, or the one with 
the highest signal-to-noise ratio. Each mean spectrum contains spectra of N SNe, where N > 3, except the ones at 
phases of t = -1-62 and -1-64 days. Thus, for SNe Ic-bl spectra with phases between t = -1-61 and -1-65 days, we use the 
mean spectrum with the nearest phase — either the one at phase t = -1-60 days or the one at phase t = -1-66 days — 
to replace mean spectra at phases 62 and 64. Since SN Ic-bl spectra have bro ader features and higher velocities than 
those of SNe Ic, we convolve the template with a Gaussian kernel (see Section 3.21 and doppler blueshift it using the 
relativistic Doppler formula. 

There are four parameters in our fit: a parameter for the width of the convolution Gaussian (ct), a velocity parameter 
(u) for the doppler shift, an amplitude parameter (a), and a wavelength-range parameter (Aw), by which we add to 
or subtract, times two, from the initial guess for the wavelength range of the SNe Ic-bl spectrum over which we will 
fit the template. The initial guess of the wavelength range is based on finding the local maxima on the blue and 
the red side of the feature and using the range in-between. We choose to use a single parameter instead of two 
(i.e., one for each end) to control the range because we want to limit the dimensionality of our parameter space, as 
a higher dimensional parameter space induces higher uncertainty in the parameter distributions, and requires more 
computational power. We use the following priors on our parameters: uniform prior for absorption velocity with 
respect to the SN Ic template (0—20,000km s“^, where positive values refer to blueshifts, i.e., negative velocities ), cr 
(0—20,000 km s“^, where positive values refer to wider kernels), and amplitude (0—3). We used a gaussian prior for 


the wavelength-range parameter Aw, where we centered the Gaussian on zero and used 100 A for the 99.7th percentile 
region. We note that these choice are only valid to up to tymax < 60 days, as afterwards the SNe Ic-bl spectra become 
too similar to the SN Ic template spectra. 

In summary, our likelihood function is: 


wri—Aw 

\np{ficbi\fic,errficu,cr,v,a,Aw) = -0.5x ^ 

x—wbi-\-Aw K 


j {ficbijx) - [ax D{v) X {G{s) * fic)]{x)}'^ 

errficbi{xY x (iV - 4) 


■ ln[err//cbi(a:)^ 


(Al) 


where ficbi{x) is the flux of a SN Ic-bl spectrum at wavelength x, D{v) is the relativistic doppler shift at velocity v, 
G{s) is the Gaussian function with width cr, fjc is the flux of the SN Ic mean spectrum, errjjcbi is the uncertainty 
of ficbh N is the number of data points between wavelength wbi + Aw and wri — Aw, and 4 repre sents the number 


of paramet ers in our model. We sample the posterior distribution function via the emcee package (Foreman-Mackey 
et al. 2013[) which implements an affine-invariant ensemble Markov chain Monte Carlo (MCMC) sampler (Goodman 


& Weare 2010). The velocity parameter v and its error bar are based on its marginalized distribution: the Doppler 
velocity parameter v corresponds to the 50th percentile in the marginalized v distribution (i.e., median), while its 




























26 



Fig. 11.— One- and two-dimension al p rojections of the posterior MCMC prob ability distributions for the spectral fit parameters for 
PTF lOqts at tYmax=0, shown in Fig. |10[ This figure was made with corner.py < |Foreman-Mackey|2016| l. 


error bars correspond to the 50*^ — 16*^ and 84*^ — 50*^ percentiles of the marginalized distribution of v (Fig. 11). We 


compute the final Fe II A5169 velocity for SNe Ic-bl by adding the Fe II A5I69 absorption velocity of the SN Ic template 
to the velocity parameter value from the fit. We compute the final uncertainty in the Doppler velocity parameter v by 
adding in quadrature the MCMC-derived uncertainty to that based on the velocity error in the corresponding SN Ic 
template. The velocity error in the SN Ic template is based on the fact the template spectrum is generated by spectra 
of different SNe Ic and at slightly different phases - thus the template velocity error was determined by taking the 
standard deviation of velocities measured in the individual SN spectra that were used to generate the SN template 
spectrum. Note that a potential mismatch between the phase of the SN Ic-bl spectrum and that of the mean SN Ic 
spectrum that was used for the fit does not impact the final absorption velocity measured for the SN Ic-bl spectrum. 

In Figure [TOl we demonstrate our template-fitting method for the SN Ic-bl PTFlOqts at phase tvmax = 0 days, using 
the SN Ic mean spectrum at phase tymax = 0 days as a template. Since the SNe Ic mean spectra are constructed 
using flattened spectra, we use the flattened version of the SNe Ic-bl spectra as well. We note that for SN 1998bw, its 
blue Fe trough is stronger than the one in the SN Ic template, at all phases. The same applies to SN 2002ap, but only 
for up to tvmax = 10 days. 

We tested our template-fitting method by applying it to SNe Ic spectra and comparing the resultant Fe II 5169 
velocities with those from line identihcations. We find that the Fe II 5169 velocities from the two methods are 
consistent within their error bars. Thus, we can be confident that using the template-fitting method, which we use 
for SNe Ic-bl does not introduce a systematic bias in the velocity measurements. We note that for SN Ic 2013dk, 
the three iron lines were blended into one feature, such that we used the template-fitting method to measure its 
absorption velocity. Our code to measure absorption and width velocities for SNe Ic-bl (and for any other SNe with 
high line-blending) is available on github under https://github.com/nyusngroup/SESNspectraLib . 


B. ABSORPTION VELOCITIES OF INDIVIDUAL SNE 


While we presented the bulk properties of SNe in the above sections, readers may be interested in the velocity 
evolution of individual SNe. Thus, in Fig. we present the Fe II A5169 absorption velocities of individual SNe in 
three panels for the different SN subtypes. For the plot that shows SN-GRBs, we also distinguish between SN Ic-bl 
connected with LLGRBs vs. SNe Ic-bl connected with HLGRBs, as well as the radio-relativistic SNe 2009bb and 
2012ap. 
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Fig. 12.— Fe II A5169 absorption velocities for individual SNe, grouped according to SN type. For the top plot that shows SN-GRBs, we 
also distinguish between SN Ic-bl connected with LLGRBs (filled circles, SNe 1998bw, 20031w, 2006aj, 2010bh) vs. opposed to the rest of the 
SNe Ic-bl (open circles). The latter include those that are connected with HLGRBs (SNe 2003dh, 2009nz, 2013cq) and with intermediate- 
luminosity GRBs (SNe 2012bz, 2013dx) and those that did not have an accompanied GRB, but exhibited very strong radio-emission, 
suggesting a relativistic engine (SNe 2009bb, 2012ap - indicated by an extra black circle around them). 


















